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Abstract 

A land cover-vegetation map with a base classification system for remote sensing 
use in a tropical island environment was produced of the island of Hawaii for the 
State of Hawaii to evaluate whether or not useful land cover information can be 
derived from Landsat TM data. In addition, an island-wide change detection 
mosaic combining a previously created 1977 MSS land classification with the TM- 
based classification was produced. 

In order to reach the goal of transferring remote sensing technology to State of 
Hawaii personnel, a pilot project was conducted while training State of Hawaii 
personnel in remote sensing technology and classification systems. Spectral 
characteristics of young island land cover types were compared to determine if 
there are differences in vegetation types on lava, vegetation types on soils, and 
barren lava from soils, and if they can be detected remotely, based on differences 
in pigments detecting plant physiognomic type, health, stress at senescence, j^eat, 
moisture level, and biomass. Additionally, literature of mapping in Hawaii was 
reviewed since national mapping systems for remote sensing do not include 
tropical island environments. Geographic information systems (GIS) and global 
positioning systems (GPS) were used to assist in image rectification and 
classification. GIS was also used to produce large-format color output maps. An 
interactive GIS program was written to provide on-line access to scanned photos 
taken at field sites. 

The pilot project found Landsat TM to be a credible source of land cover 
information for geologically young islands, and TM data bands are effective in 
detecting spectral characteristics of different land cover types through remote 
sensing. Landsat TM is the most powerful satellite sensor available to date; none 
of the current satellite sensors are ideal for the level of land information desired 
or needed in a mountainous, tropical island environment with cloud forests and 
erupting volcanoes. The environment of Hawaii requires extensive field work 
because there are established vegetation mosaics and numerous transition zones 
which change over very short distances, unlike ecozones in the temperate zone. 
The mountainous terrain covers over 4,000 sq. mi., ranging from sea level to the 
summit of Mauna Kea at 13,796 ft. along a mere 17 mile long transact. -Land 
cover change is based on climate, elevation, moisture, with primary succession 
after volcanic eruptions. Large agriculture field patterns were resolved and 
mapped successfully from wildland vegetation, but small agriculture field patterns 
were not. Additional processing was required to work with the four TM scenes 
from two separate orbits which span three years, including El Nino and drought 
dates. Results of the project emphasized the need for further land cover and land 
use processing and research. Change in vegetation composition was noted in the 
change detection image. 



It is hoped that future satellite sensors and other remote sensing instrumentation 
will address the orbit and scene constraints by developing systems that are able to 
scan a whole island in 1-2 North-South scenes and continue to improve 
bandwidths to address the plant community level, smaller land use patterns, fire 
and acid rain from volcanic eruptions, and include radar as a band, ideally with 
the same spatial resolution, for cloudy environments. 
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PREFACE 


Since Hawaii has a growing interest in using remote sensing to monitor land cover, 
land use, and environmental questions, the summary for this pilot project is 
expanded to include a wider audience than originally planned for NASA. The need 
is great on all islands, where resource managers and researchers are without tools to 
monitor the resources they are responsible for. Resource management funds are 
limited and focus is on maintenance and control needs such as fencing, cattle guards 
to protect land from pigs, cattle, and goats, bio and chemical controls. Funds are 
needed to map change in land cover and land use. 

The funds to produce this map were minimal and the result is a map which is 
considered base in areas where field data needs processing or collection. 
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Introduction 


The NASA Mission to Planet Earth has a program for environmental monitoring 
with remote sensing technology. The Environmental Analysis and Applications 
Program awarded a federal grant to The University of Hawaii's Office of Technology 
Transfer and Economic Development (OTTED), administrator of the grant for the 
State of Hawaii Office of State Planning, (OSP), Department of Land and Natural 
Resources (DLNR), University of Hawaii Planetary Geosciences, and Department of 
Geography. A portion of the grant provided for a remote sensing pilot project using 
NASA Landsat Thematic Mapper (TM) digital data with ER-2 (U-2) photography 
from NASA at Ames Research Center. The purpose of this pilot project was to train 
State of Hawaii personnel in uses of remote sensing by helping them produce a 
baseline land cover, vegetation, and land use map of the Hamakua Coast area and 
the Island of Hawaii, with a subtropical land cover-use classification system for 
Hawaii’s island environment. Each agency can further process the data for their 
agency needs. 

The pilot project produced a reconnaissance level land cover and land use map by 
pre-processing and classifying four Landsat TM scenes. Map compositions of the 
island-wide mosaics were produced of the TM raw data, the spectral land cover 
classification, and a change detection image; which combined the Landsat TM^with 
MSS data from 1977. The remote sensing software ERDAS/Imagine v 8.1 for the HP 
was used in combination with ESRI/Arclnfo v 6.0 - 7.0 for GIS. The migration of 
remote sensing and GIS to the desktop, combined with network and Internet access, 
make this technology a flexible and affordable tool for resource management, 
research, and education. 

Problem Definition 

The Hawaiian Archipelago is the most isolated major island chain on earth. The 
islands are actually the emergent tops of huge volcanic mountains, ranging from sea 
level to 13,796 feet (G. Walker). From the town of Hilo at sea level to the summit of 
Mauna Kea at 13,796 feet, road mileage is a mere 44.2 miles of mostly paved roads, 
with ecozones ranging from tropical rain forests and woodlands to alpine deserts. 
The distance from sea level to the summit is a mere 17 miles on a NO A A chart. 
Mapping land cover and land use in this tropical and mountainous island 
environment presents problems not found on continents or flatland areas with 
little topographic variability. Operating funds for field work have been declining for 
several years and ways to extend field work are needed. Therefore, alternative 
reliable methods of monitoring land cover and land use are needed for planning, 
management, and research at the landscape level. 

The last complete set of land use and land cover maps for the State of Hawaii was 
made in 1963 from aerial photography flown in 1950. Since that time, extreme 
changes in land cover and land use have taken place. These changes include decline 
of agriculture based on sugar cane and pineapple, deforestation, grazing, noxious 
weeds and pests from overseas, hurricanes, urbanization, fire, and on-going 
volcanic eruptions on the geologically young island of Hawaii. For the nation's only 
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tropical rain forests the crisis is urgent since trees catch rainfall, protecting the 
forested watershed and its ecosystems which provide water and cool temperatures 
to help sustain Hawaii's environment. Following the loss of native forests in the 
1900’s, there was an increased awareness of their importance for irrigation and 
fluming so important to Hawaii's sugar industry. Hawaiian forests were considered 
"protection forests" as opposed to "supply forests", and water production was 
considered their most valuable product (Judd 1927; Cuddihy and Stone 1990). 

The largest sections of undisturbed native communities remaining today are 
primarily Montane Wet Ohia or Ohia-Koa forest communities in areas generally 
considered unsuitable for timber production, agriculture, or ranching (Jacobi 1990). 
There is limited time available to find and implement solutions to monitor change 
in land cover for management, research, and conservation. The survival of native 
ecosystems depends on adequate, protected habitat, and research and management 
to prevent further loss of plant communities (Gagne and Cuddihy in Wagner et al). 

Purpose and Need for Action 

Land use needs to be monitored with the on-going ecological and economic crisis in 
Hawaii. Management needs visual tools and ground field site data for monitoring 
and planning land use and zoning changes in agriculture, conservation areas, 
ranching, forest reserves and plantations, to locate indigenous community and 
archeological sites, and urban encroachment. 

Hawaii contains the highest number of native plant species for any plant region in 
the world, with 89% being endemic or found only in Hawaii. Today, Hawaii's major 
ecosystems are in crisis from the continued introduction of non-native plants from 
overseas and the destruction of natural habitat which opens up sites for 
colonization by alien plants. Non-native plants amount to 47% of Hawaii’s 
vegetation at the species level, higher than first anticipated and more is expected, 
especially as Kauai is further surveyed (Wagner et al). 

The source of rapid ecological change in 'paradise' is the islands' natural beauty and 
the ever increasing number of people bringing in non-native species, intentionally 
or accidentally within the last two hundred years, since the rediscovery of Hawaii by 
Europeans in 1778. Humans thus effectively broke the natural isolation barrier of 
the Hawaiian Islands. (Mueller-Dombois; Nelson). Non-native plants arrive in 
Hawaii via the pathway of uninformed or careless inspections and shipping, or 
smuggling: airline passengers 27%, First Class mail 24%, cargo 16%, military'! 3%, 
foreign inspection 13%, and private yachts 6%. (DLNR). 

Although Hawaii accounts for 0.2% of the nation's land area, 75% of the U.S. 
recorded extinction's are Hawaiian. Today Hawaii is home to 40% of the nation's 
endangered or protected plant and animal species. (Department of Land and Natural 
Resources). Research has shown extinction is preventable if appropriate habitat 
conservation, education, and management methods are taken (Wagner et al). 
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A pplying Remote Sensing - Ecosystem Management and Research 


Incorporating satellite and low altitude remote sensing imagery with other land 
information provides a visual tool for monitoring change in landscape. The timing 
in Hawaii's history is both critical and opportune for finding ways of monitoring the 
complex interactions of both native and disturbed environments. The unique 
visual and spectral data provided by remote sensing technology is identified as one 
tool for monitoring landscape change. Offering repeat coverage of the same 
geographic area, multiple resolution image data is combined from satellite, sensors 
on-board aircraft, or spectra data from other types of instruments. 

By absorbing and reflecting energy in the environment, the imagery contains 
information responsive to vegetation type, pigments detecting health or stress, cell 
structure remaining at senescence, biomass, and moisture levels. Imagery also 
contains information responsive to water, substrate, land cover and use, 
temperature, and atmosphere. There are infinite applications for visual display, 
land cover and land use mapping, and other types of digital mapping. 
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Background 

The State of Hawaii had previously been involved in the NASA Ames Research 
Center's Western Regional Application Program (WRAP) from the late 1970's until 
the program was discontinued in the early 1980's. While participating in the WRAP, 
State of Hawaii personnel were introduced to the use of digital data collected by 
satellite to discriminate various classes and patterns of land use and land cover. As 
an adjunct to the application of remote sensing technology. State personnel were 
also introduced to the use of the relatively new technology of geographic 
information systems to integrate mapped information with processed satellite 
images. NASA's VICAR/IBIS was used for producing a Prime Forest Overlay from 
the supervised classification of the island of Hawaii. 

When the WRAP was discontinued in the early 1980's due to Federal budget cuts, 
advances in computer technology and GIS software eventually enabled the State to 
initiate a GIS program in the late 1980's. 

Since the State's GIS was envisioned to be a multi-agency system accessible to all 
State agencies with an interest in spatial data, the Office of State Planning within the 
Office of the Governor, was designated as the lead agency for planning and 
developing the State's GIS program. The foundation of the system wa^=the 
development of a centralized database which would be accessible to all user agencies 
to avoid duplication of effort in the costly area of database development, while 
maximizing efficiency in data management in State government. 

Although the implementation of a GIS offered a wide range of applications, the 
initial focus of the database development efforts were concentrated on land use 
planning and resource management uses. In addition to the development of the 
requisite base maps, most of the data layers which were initially developed included 
environmental and natural resources-related data layers. 

In March 1993, NASA representatives contacted State of Hawaii officials to discuss 
the potential benefits of incorporating the use of remotely sensed data with the 
State's existing GIS to support resource planning and ecosystem management 
decision making. State of Hawaii personnel expressed interest in this initiative, 
since it would augment the existing GIS database and also provide an opportunity to 
evaluate the feasibility of using Landsat Thematic Mapper (TM) imagery to develop 
land use, land cover, and vegetation data sets. _ 

In its previous experiences through the WRAP, the State's Landsat Demonstration 
Project utilized Multi-Spectral Scanner (MSS) data to develop land cover/land use 
classifications of the island of Hawai'i and O'ahu. Since the State had already 
worked with MSS data with 4 spectral bands and approximately 1 acre resolution 
(79mx79m), this new initiative offered the State an opportunity to evaluate the 
viability of using TM data which has a higher resolution (28.5m x 29.5m) and with 7 
spectral bands of covering land. 
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State of Hawaii Consortium for Technology Transfer 

Cooperative alliances between public and private organizations are critical to find 
the authority and resources to address the needs of Hawaiian ecosystems, planning 
effectively for new land use and creating new jobs, recovering and sustaining forest 
lands and the remaining island ecosystems of Hawaii and the Pacific. (Tropical 
Forestry Plan 1995). Although the State of Hawaii constituencies have their own 
priorities and responsibilities, collectively they need visual tools and base maps for 
assessing and sustaining Hawaiian ecosystems at risk. Imagery provides one layer to 
monitor the changes visibly evolving in land cover and land use. 

The University of Hawaii's Office of Technology Transfer and Economic 
Development (OTTED), representing a number of State and University 
Departments, applied for a Federal grant from the National Aeronautics and Space 
Administration (NASA). OTTED served as the administrator of the grant, on behalf 
of the State of Hawaii Office of State Planning (OSP), Department of Land and 
Natural Resources (DLNR) and the University of Hawaii's Department of Planetary 
Geosciences and Department of Geography. In August 1993, the above consortium 
was awarded a grant from NASA Headquarters to in part, demonstrate the potential 
benefits of the Environmental Analysis and Applications Program being developed 
by NASA s Office of Mission to Planet Earth, while helping the State of Hawaii re- 
establish its in-State remote sensing analyses and image processing capabilities. 

Project Requirements and Goals 

By carrying out a remote sensing pilot project using Landsat TM, to produce a 
baseline land cover classification map of the island of Hawaii, the State of Hawaii 
hopes to attain the following project goals. 

• Design a system for State employees to process and retrieve remotely 
sensed image data on the State of Hawaii GIS. 

• Train State employees to produce a baseline land cover, land use, and 
vegetation classification map of Hawaii Island from Landsat TM data. 

• Integrate processed Landsat image data with the State of Hawaii's 
existing GIS database design. 

• Evaluate the usefulness of Landsat imagery for carrying out land 

cover /land use and vegetation classifications. -»• 

• Evaluate the usefulness of Landsat data for detecting changes in land 
cover /land use, vegetation and the environment. 

• Establish a new base line of data to be used for "change detection." 

• Establish within the State of Hawaii, the capability of reformatting and 
distributing remote sensing digital data. 

• Position the State of Hawaii to become a possible "Center of Excellence" 
for the Pacific Basin in GIS and remote sensing technology. 
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Implementation of the Plan 

A collaborative effort by the State and University of Hawaii entities was initiated in 
order to meet the stated project goals. OTTED, as the grantee, was in charge of 
overall project management and federal grant administration. OTTED also served as 
the State's point of contact with NASA's Earth Observing System (EOS) Program. 

The University of Hawaii's Planetary Geosciences Department was in charge of 
receiving the acquired data and reformatting multiple sets onto CD-ROM. They also 
were the point of contact for distributing copies of the CD-ROM to State of Hawaii. 

The University of Hawaii Geography Department had students participate in the 
initial ERDAS software training and later on in the change detection portion of the 
project. Remote Sensing is offered in addition to their cartography and GIS 
curriculum. 

The Division of Forestry and Wildlife and Office of State Planning assumed the role 
of the application project team in order to develop in-house image processing 
expertise by actually participating in the image analysis and classification pilot 
project for the Island of Hawaii. The Division of Forestry and Wildlife also provided 
ground support for field data collection. _ 

OSP served as the State agency liaison with OTTED and was also responsible for 
upgrading existing State GIS equipment and procuring image processing software 
for State agency use. DOFAW was in charge of coordinating and leading the 
necessary field work efforts on the Island of Hawaii. 

Observations 


This technology transfer pilot project outlines methodology to incorporate and 
process remote sensing imagery in a centralized location with access to State 
databases for State and University personnel. The sensors included are NASA’s 
Landsat TM and MSS, and color infrared photography (CIR) from the ER-2 (U-2). 
The need is great on all of the islands where resource managers and researchers are 
left without access to these visual tools for management and research of the lands 
they are responsible for. It is time to expand the resources and knowledge to the 
locations that need them most. The interest level in remote sensing is seeded, as 
exemplified by the recent purchase of statewide coverage of SPOT Image data. Public 
and private personnel will form alliances to use and share data and information. 

Acquisition of Hardware and Software 

A portion of the grant was allocated for the purchase of workstation hardware and 
software designed to process remotely sensed image data and to integrate the data 
into the state's GIS. Remote sensing software provides tools to process image data 
with an overview of landscape (30m), where land features look small; or, very high 
resolution data (l-10m) from an airborne scanner, where land features look big. A 
total of 2 ERDAS/Imagine software licenses were purchased for the State GIS. 
ERDAS is the industry leader in remote sensing software and provides the most 
compatibility with the State’s Arc/Info GIS. 
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Hiring of Consultants 

Since the State and the University of Hawaii did not have personnel available to 
produce the pilot project, and in order to effectively carry out the technical 
requirements of the grant, the State opted to use a portion of the federal grant to hire 
a consultant to facilitate the technology transfer process. In this regard, OTTED 
solicited proposals to develop in-State expertise in analyzing and processing 
remotely sensed data and then applying such information to the State's GIS. 

A proposal submitted by Geographic Decision Systems International (GDSI) and 
Hogan Co. was selected. The proposal offered project management at the State site 
using the systems approach to a remote sensing project. This approach guides image 
database organization and parameters for pre-processing data for the sensors used in 
the project. The remote sensing project system is composed of three functional 
dimensions which correspond to the ecosystems of the project site. The remote 
sensing project design is further divided into subsystems for processing: input, 
analysis, and product database subsystems. They are represented by the project flow 
on the following pages. 

The remote sensing phase of the project was led by Hogan Co., and the GIS phase 
was led by GDSI. GDSI provided custom software where needed for map production 
and the multi-media AML to view photos of GPS sites. Together with Office of State 
Planning and DLNR they provided on-going support for map production, GPS 
editing, and the multi-media AML program LSVIEW. 

Training in the basics of image processing and GIS were taught using the State's 
hardware/ software and the pilot project's data, to produce a baseline land cover, 
land use, and vegetation map. Small scale maps used for administrative purpose 
and large scale maps for the field were produced. This approach delivers the 
capability of further processing, which allows State personnel to process the data in a 
working environment they are familiar with. Independence is established to 
analyze remote sensing imagery, produce maps, and to establish remote sensing data 
within the GIS database design. Introductory training of 17 participants included 
OSP, DLNR, ICSD, and UH. Following the initial two weeks of training, the 
consultants worked closely with the State personnel when they were available to 
participate in the project. 
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LANDSAT TM - PILOT PROTECT SYSTEM DESIGN 

The project flow was as follows: 

• Introduction by NASA 

• A Remote Sensing System - scene, sensor system, processing system 

• State of Hawaii Requirements and Objectives 

• Establish feasibility 

• Project Plan 

• Project Implementation 

• Remote Sensing Classification hierarchy -major land cover ecosystems 

• Output Database Subsystem - processing standards for Landsat TM in 
State GIS, archive /retrieval format/media, geometric grid: terrain, cell 
size for image/GIS layers, map scales, interactive use, classification tree. 

• Input Subsystem - scene acquisition and formatting for distribution 

• Deposit new TM and 1977 MSS raw data on GIS file server _ 

• Import data into ERDAS/ Imagine for image modeling and processing 

• Analysis Subsystem - georectify TM data to UTM Zone 4 and correct for 
terrain displacement with 32m pixels, to comply with the State GIS 

• Unsupervised classification of each of the four scenes for Hawaii Island 

• Ground field data collection 

• Produce GPS Point Coverage of ground field data 

• Supervised classification and training field extension of north scenes 

• Field data collection and verification 

• Supervised classification with training field extension of four scenes 

• Produce whole-island raw data and classification mosaics by modeling 

• Unsupervised classification examples on whole-island raw data mosaic 

• Change Detection Analysis - georectify 1977 MSS to TM in State GIS 

• Change Detection - of combined MSS_TM image 

• Supervised Classification Masking Examples in ERDAS /Imagine v8.2 

• Vectorize supervised classification for GIS layer in Arc/Info Version 7.0 

• Map production / GIS support as maps produced for Output Database 

• Create multi-media application LSVIEW to demonstrate use of GPS 

• Products archived on: tape, optical disk, film, CD, paper, or file server 
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Technology Transfer Output Product Design 

Building an image database requires agreeing to the same standards for product 
design, reassessing them as needed, and documenting sources on image/GIS layers. 
Specific needs must be addressed in the process if land cover information is 
incorporated into a larger mapping system of resource data. 

Map products for the pilot project were produced with ERDAS/Imagine and 
ESRI/Arclnfo software, but Arclnfo was used more by the State since they have 
custom AML programs in place to automate map production. The University of 
Hawaii used the SUN platform with demo license of ERDAS/Imagine v8.2 to 
produce their products. The products from this project are outlined below and listed 
in the Appendices as Appendix H: 


1) an unsupervised land cover classification of the Hamakua test site area, 
extending from west of Honokaa to Laupahoehoe. 

2) a base subtropical land cover-use classification system for use with remote 
sensing data in Hawaii's island environment. 

3) a base-land cover, land use, and vegetation map of the island of Hawaii^ 
produced at 1:250,000 in image, hard copy, and 4X5 negative format. 

4) a digital six band color mosaic of the four TM scenes including Maui. 

5) color composite of TM bands 7-5-3 scaled at 1:250,000 showing separation of 
vegetation from re vegetated lava flows, lava, agricultural areas and urban 
classes (Hilo, Kona, Waimea) where cloud cover permitted. Digital, hard 
copy, and a 4X5 negative transparency is available. 

6) a ten band change detection composite image of Hawaii with multiple 
resolution Landsat data combined from 1977 MSS channels 4-5-6-7 as layers 
1-2-3-4 of the image and TM bands 1-2-3-4-5-7 as layers 5-6-7-8-9-10 of the 
image. Digital and hard copy examples of areas showing significant change 
are in 8.5 X 11.0 inch format and a digital composite of TM 7-5 with MSS 6 
as layers 10-9-3. 

7) Vectorized land cover, land use, and vegetation maps formatted in 7.5' 
quadrangles scaled at 1:24,000. ( 1.5gb data can be loaded as needed). 

8) State personnel with training in remote sensing project system design. 


Both electronic and hard copy maps are available, located at the State GIS for 
manipulation and further processing by the agency personnel. Image, vector, CD, 
hard copy maps, and 35mm and 4 X 5 in. negative transparencies were produced. 
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Cost Effectiveness of Remote Sensing Surveys 

This is a value judgment area that needs reassessment throughout a project as 
factors effecting project flow are identified. If the level of information awareness or 
need increases, the cost and value of the map also increase. Readjusting cost and 
value of the map is required. 

Vegetation patterns mapped in an area with minimal field samples, levels of data, 
or funds is unfortunately reflected on the map. With powerful sensors like Landsat 
TM and use of multiple resolutions and scales of data, techniques for data mining' 
are available to discern multiple levels of land cover information, depending on 
need. 

The Island of Hawaii - Project Site 

Geologically, Hawaii is the youngest island in the archipelago and expanding with 
the addition of new land from the active lava flows of Kilauea's on-going eruption 
1983-present. The big island is (4,638.2 sq.mi.), two times the combined size of all 
other Hawaiian islands, four times larger than the eastern state of Rhode Island, and 
approximately the size of Connecticut (4,845 sq.mi.) with Kilauea forming new land. 
Hawaii has well-defined vegetation zones based on moisture-climate, elevation, and 
vegetation-habitat zones with physiognomic types of, e.g., forest, shrub, grass. (Jacobi 
et al. 1986; Gagne, and Cuddihy ). 

Hawaii is the largest island in the Hawaiian Archipelago. The archipelago is formed 
by volcanic islands consisting of 8 major islands, islets, reefs, and small atolls in the 
central Pacific. The archipelago extends from Hawaii Island in the southeast to Kure 
Atoll and Midway in the northwest; totaling 10,932 sq. mi., with a land area of 
6,459sq.mi. (16,729 sq.km.). The 8 main islands comprise more than 99% of the land 
area and extend from Hawaii Island at 18 degrees N latitude in the southeast to 
Niihau at 22 degrees N latitude in the northwest, and the northwest islands and 
atolls extend to 28 degrees N latitude (G. Walker). 

Mapping land cover and land use in this tropical and mountainous island terrain 
presents problems not found in continental, flatland areas with little topographic 
variability. The uniqueness of the vegetation is a result of geologic history, extreme 
isolation from continents, variation in substrate, topography, and climatic changes 
over short distances, which allow the vegetation to migrate radically into available 
niches. These factors contribute to the diverse mosaic formations of the natural 
plant communities. (Wagner 1994, Mueller-Dombois, Bridges, & Carson, 1981). 

Physical Ecosystems of Hawaii Island 

Hawaii Island is formed by 5 volcanoes: Kohala in north, Hualalai in west, Mauna 
Kea (13,796 ft. or 4,205m) and Mauna Loa near the center, and Kilauea on the 
southeast slope of Mauna Loa. Kilauea is the most active volcano in the world with 
the current eruption expanding 1983 to present. Kilauea is now the largest and 
longest rift eruption in recorded history, forming Pu'u O'o on January 3, 1983. In 
1986, a change in eruption activity formed the shield Ku'paianaha (the mysterious) 
with a lava lake within it's summit feeding a network of lava tubes, which flowed 
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to the ocean producing 500,000 cu. yds. of molten lava daily to 1992. A fissure opened 
in Pu'u O'o and this fissure along with others continues to feed lava tubes 
extending to the sea from 1992-date (HVNP). 

Mauna Loa (13,677 ft. or 4,169m) is the world's largest volcano, rising almost 30,000k 
from ocean floor and approximately 70mi. (97 km) wide at its base, erupting recently 
in 1950, 1975, with the last eruption in 1984, with subsequent hot spots visible on 
infrared imagery to date. 

The last eruption of Hualalai was in 1800-1 and the last eruption for Kohala was 
some 60,000 years ago. 

There are diverse climatic and precipitation changes over short distances from 
topography ranging from sea level to the summit of Mauna Kea at 13,798 ft. 
elevation. This elevation change is a mere 17 mile long transact when drawn on a 
NOAA navigation chart of the Hawaiian Islands (ed. 1991) or 44.2 miles by roads. 
Recent geologic age, solar radiation and latitude, slope, and aspect affect veg type and 
density which causes change in spectral values. 

Biological Ecosystems of Hawaii Island 

Vegetation is the obvious biological process visible in the imagery. The Hawaiian 
flora is well known for being the most distinct plant region in the world, dominated 
by native plants found no where else. The native land cover of Hawaii are forests 
and Hawaii is the only state with tropical rain forests and habitat (Cuddihy and 
Gagne). 

Island ecosystems are different than continental ecosystems due to geographic 
isolation, predominately small geographic areas, and recent geologic age with well 
developed biological differences resulting in extremely visible interaction between 
native species and native with non-native species (D. Mueller-Dombois). Forest 
succession is on geologically young land, there are established vegetation mosaics, 
and abundant transition zones with diverse vegetation types changing in short 
distances and time spans, based on climate, elevation, and moisture. 

Gagne and Cuddihy modified Jacobi’s outline of vegetation zones, based on five 
elevation zones of coastal, lowland, montane, sub alpine, and alpine. Three 
moisture zones are described as wet, medic, and dry, with five physiognomic plant 
types defined as grassland, shrubland, forest, open forest, and parkland (Gagne, and 
Cuddihy). The combined information was adapted to fit the NASA pilot project 
based on the data's spatial, spectral, and temporal resolution (including two orbits 
and three years) and includes Tropical Wildland, Physical, and Manmade systems. 

The diverse ecosystems are complex based on elevation, moisture, and climatic 
zones, volcanic activity, and time described here for Landsat TM’s resolution and 
the field work collected. The minimum resolution is 6(32. 0m x 32.0m) spatially, 
derived from the six Landsat spectral bands, over three years time including El Nino 
and drought years, with atmospheric effects of clouds, precipitation, ice, snow, and 
VOG. 
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Manmade Systems - Land Use and Vegetation Introductions 

Hawaii has abundant indigenous populations and ancient archeology sites which 
interact in both remote and urban environments. Colonization of Hawaii by the 
Polynesians is recorded from archaeological sites with settlement time 300-400 A.D. 
at Bellows Waimanalo on O'ahu (Kirch 1974; Cuddihy and Stone 1990). By the 6th 
century, all islands had permanent settlements. Alien plant introductions began 
with the Polynesians, who brought with them about 32 species. Of these 32 species 
less than 25 escaped cultivation (Cuddihy 1990; Nagata 1985)). The enormous 
increase in non-native plant and animal introductions has occurred in the past two 
hundred years by humans. (D. Mueller-Dombois). 

Economic changes, like ecological changes, are more apparent on islands than 
continents, again because of isolation and limited resources on generally small and 
geologically young land masses. 

There are major economic changes present with the immediate end of the sugar 
industry and 80,000 acres of abandoned cane fields in transition. The last harvest for 
Hamakua Sugar taking place during the project's field work in 1994 and in the 
winter of 1996 for Kau, C & H Sugar in the Pahala area. The Kahala Sugar Company 
discontinued production in 1977, the same year of the Landsat MSS data. Some4and 
use conversion is on-going. There is a potential for reforestation with native Koa, 
non-native forest species of eucalyptus, sugi, ash, or bamboo, or to extend crops such 
as truck farming, wine grapes, coffee, taro, macadamia nut, cattle, and nurseries of 
flowers and trees. 

Island of Hawaii - Remote Sensing Project System 

A remote sensing project system is composed of image scenes for a project site, the 
resolution domain of a sensor system to record information, and a processing 
system for analysis and output products. This system corresponds to physical, 
biological, and manmade ecosystem domains of the project site. The scene 
corresponds to ecosystems in front of the sensor including the atmosphere, the 
multispectral sensor records visible, infra-red, and thermal information contained 
in the scene, and the processing system analyzes the information visually and 
numerically, and produces trained personnel and image or map products (See 
sections on: Landsat TM - Pilot Project System Design and Product Design). 

For example, irt this project image data from three dates of one sensor (LandsaLTM), 
is combined with ancillary information to produce the multidate classification map. 
Image data from more than one sensor is combined to produce the multiple 
resolution change detection composite mosaic of Hawaii, with Landsat MSS and 
TM. 
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Selecting Landsat TM Scenes for the Island of Hawaii 

In addition to the federal grant, NASA provided seven Landsat TM scenes as part of 
their bulk purchase plan of data. The scenes of Kauai, Oahu, and Hawaii were 
unacceptable due to clouds over the island land masses and could not be used for 
land cover classification or change detection. The data was ordered in the Space 
Oblique Mercator (SOM) projection with a cubic convolution resampling method. 

This project typifies image processing problems found in tropical island 
environments at a latitude where the sun angle and azimuth are close. Imagery 
from different dates was used because of the separation of adjacent satellite paths (N- 
S orbit) and cloud cover on the available data. Visible atmospheric conditions occur 
over cloud forests and from volcanic activity (VOG). Forest succession is on 
geologically young land. 

For the pilot project site of the island of Hawaii, four TM scenes were required to 
mosaic a whole island composite image. A planning meeting was held in May 1994, 
to meet the project team and select big island data from 8.5 X 11 in. black and white 
copy. The quality of available imagery was an immediate concern for vegetation 
classification. The scenes with the least amount of cloud cover and temporal 
differences were selected by University of Hawaii Department of Planetary 
Geosciences/Hawaii Institute of Geophysics (PG/HIG), Department of Geography, 
Office of State Planning, Hogan Co., and GDSI. 


Landsat Path / Row Index 



Figure 1 : Dates of digital scenes with corresponding Orbit Path and Row are 
indexed along with their map extent. 
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The data cover three years and include both El Nino and drought years, as well as 
seasonal variation. Mauna Kea received 17 feet of snow in March 1990, with 
snowfall recorded through July 1990. On the 1991 drought image named Kalapana, 
the visible Kilauea volcano eruption sites are from Pu'u O'o (formed in 1983) and 
the shield Kupaianaha (formed in 1986). 'Skylights' of molten lava, burning 
vegetation, and VOG are visible. Kupaianaha was active from 1986-1992 when 
activity moved back up the east rift to Pu'u O'o, where it remains active to date. The 
effects of the visible atmospheric differences and moisture stress in the vegetation 
on the drought years adds dimensions of data to process in an environment that is 
already complex. The impact of these differences on the classification are impossible 
to determine definitively (Matt McGranaghan). 

APPENDIX A - Landsat TM Scenes and NOAA data 


Sensor Characteristics - NASA's Landsat TM Sensor 


The sensor characterizes the scene by absorbing and reflecting radiation energy from 
the ground site. Remote sensing image data is acquired by a scanner with detectors 
for recording defined bands of reflective and emitted light in the visible, infra.-red, 
and thermal wavelengths. Through a process known as digitization, the analog 
signal is converted into a numerical format of binary digital data. Through the 
telemetry process, the resultant image data can be processed on a computer with the 
individual radiation bands displayed as intensity measurements in a gray scale with 
digital numeric values (DN) ranging from 0-255. Scanners record in either pre- 
defined bands as Landsat does, or they record in bands that are selectable and defined 
by a specific need to detect spectral land cover information. 

The bands are designed to detect health (i.e. chlorophyll, carotene pigments), 
senescence/ stress (i.e. carotene, xanthophyll, and anthocyan pigments), biomass 
(dense or open), moisture or dryness in vegetation and soils. Physical systems 
include water, ice, snow, and atmospheric conditions such as clouds and VOG. 
Revegetation on flows, soil types, and characteristics of lava, such as type and age are 
also detected. 

NASA's Landsat TM scans bands of electro-magnetic energy in seven wavelength 
areas, including pre-selected bands in the visible blue, green, and red; reflective^near 
and middle infrared; and the thermal infrared regions. Various band combinations 
are used to discern the complexity of azonal substrate regions and zones of 
vegetation-substrate cover observed. TM bands are used for analyzing the separation 
of lava, soil, vegetation moisture-heat, and vegetation pigmentation concentrations 
characterizing health or stress. These bands also detect the heat from the skylights' 
of active lava flows and older, cooled lava, atmospheric effects, and biomass. 

Landsat TM data is digitized to 8-bit precision per pixel and is scaled with 256 gray 
levels of radiometric data displayed in a black-and-white gray scale ranging from 
black = 0 and white = 255. The data were ordered in Space Oblique Mercator map 
projection, with earth ellipsoid of International 1909. The data was resampled by the 
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cubic convolution method, producing a spatial resolution size of 28.5m x 28.5m per 
pixel. The format was band sequential or BSQ, with the data of each TM spectral 
band organized in a separate file. 

The Pilot Project Processing System 

The processing system for this project is located completely in Hawaii and includes 
the main processing areas of georectification, classification, mosaicing, change 
detection, and map production. The scenes were georectified to UTM Zone 4 to 
comply with the State database design. Land cover, vegetation, and land use 
classifications were produced. A change detection image was produced by combining 
the Landsat TM and 1977 Landsat MSS. Map products were produced with 
ERDAS/Imagine and ESRI/Arclnfo software and are described in the Technology 
Transfer Product Design Section, Report on GIS Support, and in the Appendices. 

Preliminary Processing of the TM Scenes 

The University of Hawaii Department of Planetary Geosciences/Hawaii Institute of 
Geophysics (PG/HIG) was in charge of pre-processing the raw digital TM data files by 
formatting them for distribution. The data were delivered to PG/HIG from EROS on 
high density 9-track tapes in band sequential (BSQ) format. Hogan Co. workecLwith 
PG/HIG to make sure the data entry process would be as simple as possible, since 
this would occur in a secured room without access to project participants. PG/HIG 
stored each image band in a separate file, with the associated header and calibration 
information stored in separate text files. The reformatted data were provided to 
agency participants on high-density 9-track tape and CD-ROM. 

Loading Image Data to the State GIS File Server 

Training began in July, 1994 following the June arrival of ERDAS/Imagine v 8.1 for 
the HP platform, which the State OSP has. An overview and history of remote 
sensing was presented for 20 State agency personnel from Office of State Planning, 
Department of Land and Natural Resources - Division of Forest and Wildlife, and 
Information and Communication Systems Division (ICSD), where the State GIS lab 
is located. The first week of initial training focused on the basics of remote sensing 
with students using the images provided with ERDAS/Imagine software in the first 
few lab session until the importing of the project data was possible. 

Ideally, raw data is loaded directly into Imagine from EOSAT or EROS source 9^track 
tapes or CD-ROM using the import menu function. Because the initial release of 
Imagine v 8.1 lacked peripheral support on the HP-UNIX platform, the import 
function did not work with the State's hardware. 

First, the project team attempted to read the EROS 9-track tapes directly into Imagine 
using Imagine tape import function. This was not recommended for the duration of 
this project because the tape drive was in a secured environment and it is 
unsupported by ERDAS. Several phone calls to ERDAS hardware support in Atlanta 
enabled reading of the first tape volume of three. A bug in Imagine BSQ import did 
not prompt the system manager to mount the second and third tape volumes of the 
multivolume NW big island scene. 
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Data entry resorted to CD-ROM in the UNIX operating system command tar', 
which was a distribution method suggested as a test in the proposal for this project. 
This turned out to be the only way data could be easily transferred from the 
University of Hawaii PG/HIG to the pilot project GIS lab site. 

Each image band and calibration file were archived as separate files on CD-ROM and 
read to the file server using the ’tar' command. The Imagine Viewer does not 
recognize these files since they do not have a ERDAS Imagine image file extension 
(image band.img), and they cannot be viewed as a single gray scale band or as a color 
composite on a monitor until they have this extension. 

Each band was read into Imagine using the import function with the generic binary 
read option, to obtain an Imagine image file extension. This is not an easy method 
of handling image data for newcomers because the data is viewed in numerical 
format, and it can be difficult to see where a leading header record and trailer end 
and the image begins. The text files contain the image dimensions which is needed 
to import data. Initial disk space allocation was insufficient to load the NW scene for 
several days, but in a few days room was made for the two north scenes of the 
island. Individual gray scale bands of Landsat TM data could be viewed on the 
monitor, with the 256 levels of spectral data ranging in a scale of values from 0-255, 
black equal to Q and white equal to 255. Values in between are viewed in shades of 
gray. A color composite consisting of all the imported bands still needed to be made 
for processing the image data. The viewer menu and button options did not have a 
way to build a composite image from the individual gray scale bands and a spatial 
model was needed to generate a color composite consisting of all six bands. 

Creation of Composite Imag es 

The Spatial Modler provided the graphics tools to build custom models consisting of 
a raster layer for each image band, a function tool for data generation to stack the 
layers, and a color composite raster layer. Starting with the north west scene, 6- 
dimensional color composite images were produced of TM bands 1,2,3,4,5 and 7, 
using a custom stack model in Imagine. Bands 1-5, and 7 have different spectral 
bands of information, but spatially they are all 28.5m x 28.5m. Three spectral bands 
of the composite can be displayed at one time through the red, green, and blue filters 
of a color monitor. Composite processing time required 4-6 hours per scene. TM 
band 6 was made into an Imagine image file and saved on line because Band 6, a 
thermal infrared covers spectral bandwidth (10.4-12.4um) and has a different spatial 
resolution of 120m. Allocation of space consisting of two file servers became 
available over several months time, and these steps were repeated for the two 
southern scenes of the big island. 

Remote Sensing Classification Systems 

Building an image database requires agreeing to the same standards throughout a 
project and when there is a need to reassess. Knowledge of the domain of the 
sensor(s), the image data, the level of land cover information class needed, and 
nomenclature for a base classification system are required before the classification 
process. 


16 



Standards for preprocessing, georeferencing, and minimum map unit are needed. 
Image data has a specific resolution and is typically without a scale, except for 
U.S.G.S. digital orthophoto quads (DOQ) which are not available for Hawaii to date. 
Knowing the source information of data is critical in Hawaii, when data from 
different zones, island projections, and digital DLGs are combined with data from 
paper topographic maps or orthophoto quads. 

Ancillary photography for photo interpretation and signature development is used 
with vegetation descriptions and keys. Field sampling includes GPS for signature 
development during the classification and analysis process, and for field 
verification. Nomenclature is developed for the ecozone domains and the level of 
vegetation, substrate, and other land cover needed. This changes during analysis, 
the map verification process, and accuracy assessment. 

National resource systems including imagery are being developed which rely on 
field data and interpretation with the local resource agency personnel responsible 
for managing the local area. 

In general, national mapping systems for remote sensing reflect the continent, 
excluding vegetation differences of subtropical and tropical zones, especially island 
environments such as the State of Hawaii, U.S. territories, and other Pacific islands. 
Some of the vegetation and image processing differences are explained in this 
report. 

Ecosystems Visible in the TM Imagery 

Pattern recognition recognizes multiple spatial and temporal scales in ecosystem 
processes (Allen and Starr 1982). Ecozones in spectral imagery range from general to 
specific, depending on resolution of the data (spatial, spectral, temporal, and 
atmospheric), information class need, or by combining multiple resolutions and 
scales of digital imagery with ancillary photography, GPS, and other field site data. 

During training, the land cover of the project site was divided into two broad 
categories based on cloud free land cover seen in the imagery and needed on the 
map (signal), and data to generalize because it is not visible, needed, or accounted for 
(noise). An example is water, some agency participants wanted water on the map 
and some did not. Water classes were obtained from an automated classifier and 
saved as a file which can be used to mask water for subsequent classifications. This 
approach is used in multi-agency projects with different priorities and with project 
sites containing complex systems such as Hawaii. 

After dividing the scene into the two groups, the visible land cover was further 
divided for the initial classification scheme. Three top level subtropical ecosystems 
were listed: Physical, Biological, and Man-Made Systems. Subsystems reflecting the 
30m TM resolution and the sub-tropical island environment were listed under each 
top level system, along with the legend from a coarse level classification scheme for 
1977 Landsat MSS data with 57 x 79m resolution. The classification system that 
evolved is adapted for using the Landsat TM sensor in Hawaii's island 
environment, and based on the processed GPS and non-GPS field data. 
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Island of Hawaii - Remote Sensing Classification System 

The classification system that evolved in the pilot project is for land cover resolved 
from Landsat TM remote sensing data, collected field data, and resolved adaptations 
from Hawaii vegetation classification systems. Vegetation from the classification 
system of Gagne and Cuddihy, and Jacobi Level II, were combined with 
modifications to resource oriented remote sensing systems to account for the 
differences in a subtropical island environment. 

Remote sensing systems evolved from the first land use system for remotely sensed 
data developed by Anderson (U.S.G.S., 1974), with consideration of a nomenclature 
framework for different agency use and ecosystems from (Rhind and Hudson, 1980), 
and consideration of tropics systems and islands from life zone ecology for Costa 
Rica (Holdridge 1966,79). The relationship between level of information class detail 
in the classification system and sensor systems is from (NASA, 1983 Botkin et al., 
1984). 

The development of the remote sensing classification system was based on the TM 
image resolutions, the collection of field site data, and the software. The minimum 
image map unit for area measurements and location control is based on the pixel 
resolution domain of Landsat TM 6(32. 0m x 32.0m) grid cells, which means^and 
covers smaller than this area are not classified. A minimum image map unit of 10 
to 40+ acres is more common and valid. Exceptions to this general rule that are 
often cited for the purpose of example are streams, gulch and riparian vegetation, 
pixels for monitoring change and transition, roads, or bridges (i.e. Hilo breakwater 
and Saddle Road). 

For the pilot project, a hierarchy of desired land cover information classes was 
constructed in information-tree or taxonomic form, descending from general to 
more specific classes. As field data was acquired during the project, classification of 
the hierarchical tree was expanded in descending order. 

Top Level Ecosystems 

These systems are presented in a hierarchical information tree form and subdivided 
by criteria from the three main ecozones processes seen in the project site: Physical, 
Biological, and Man-Made Ecosystems. These systems are further divided by 
information class need, physiognomic plant type (forest, shrub, grass), moisture, 
elevation, cliirfate, separable spectrally or by masking, and pilot project Tunds 
constraints. Under Biological Systems, vegetation is the subsystem needed for the 
map. Native vegetation is tropical, whereas the origin of non-native vegetation is 
world-wide. Under Manmade Systems, land use plantations of sugar cane and 
macadamia nut were visible in the Lowland Zones where large fields occur. Nursery 
roof tops were also visible. Physical Systems are largely azonal occurring at multiple 
elevation zones, except for water classes. 

An example of change information extracted from the system is remnant forest and 
revegetated lava which were listed in the classification scheme of the last project 
and again on the first day of project participant training. Originally there were three 
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GPS points for revegetated lava, but lack of field time to collect points necessary to 
account for all the variables caused misclassification of some of these areas as cloud 
(spectral space must go somewhere). The information class Revegetated Lava was 
considered an intermediate class, needing additional ground truth to discriminate 
and label the actual information classes. 

On the second field survey, a transact was made across accessible lava flows in the 
overlap area of the north scenes to account for the enormous variability of pioneer 
land cover on the flows. Revegetated lava became several classes further down the 
hierarchical tree including the 1984 A'a flow with licken which was inaccessible on 
a previous survey. This process was repeated for accessible areas on the Puna and 
Kalapana scenes, although A'a-Lichen is colored differently in these areas. 

Some Ohia classes on the mosaic are considered intermediate classes in the Lowland 
and Montane elevation zones where mixed classes for co-dominant cover types are 
missing. Additional data, funds, and time are needed to develop Ohia - Strawberry 
Guava and Ohia - Myrica faya information classes. Also, Myrica faya is mapped with 
only the core area of concentration and a lowland occurrence with Molasses Grass. 

In many ways, this organization of classification layers represents the temporal 
differences in the data and the order they were processed in, with the least amount 
of processed field data following the elevation areas with cloud cover (Kona, Puna, 
and parts of the southern scenes). 

Processing for the temporal differences provided some control over the temporal 
differences and allowed the classification system to be conditioned on a previous 
classification, so the categories can be changed at any branch node (i.e. add 
vegetation information based on missing variability of type, elevation, slope, or 
aspect). 

Figure 2 on the following page (and Appendix C - Classification System- Charts of 
the Vegetation and Other Landcover Types in the Landsat TM Classification System, 
Signature file Class Names list with corresponding Legend used for the pilot project 
of the Island of Hawai'i). 
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Georectification 


Building a shared image database requires agreeing to the same standards for 
georeferencing or documenting sources for shared data requiring further processing. 

Geometric corrections provide a system to represent an image on a plane surface or 
grid with known map coordinates, or to align an image with another image as is 
used in medical imaging and remote sensing change detection analysis. 

Remotely sensed data are imaged by a scanner on a platform such as a satellite or 
aircraft and have distortions caused by the rotation and the curvature of the earth, 
atmospheric conditions, or by the sensor being used. The scanned image data is 
grided, and each grid cell is called a pixel with column/ row (x,y) coordinates instead 
of ground coordinates like a map. Georectification is the process of coding pixels in 
an image to the geographical coordinates of a map projection or another coordinate 
system using a nth order polynomial equation. The pixels of the new grid do not 
align or fit at corners with the pixels of the original grid, requiring a process known 
as resampling. Resampling extrapolates data values for pixels for the pixels on the 
new grid from the values of the original (source) pixels, although there are methods 
to shift and rotate pixels without interpolating new pixel values in the new corner 
areas of neighbor pixels. 

There are several considerations for the method used to georectify the projects data, 
including: map scale requirements for field work (1:24,000), the use of Global 
Positioning System (GPS) and DLG data for ground control points, map products 
(including vegetation classification and GPS overlay), the resolution of the TM data, 
existing and needed map production tools, and integration to the State GIS database 
design. (Radiometric errors of TM or SPOT satellite data are corrected by the 
companies responsible for their distribution, namely EOSAT and SPOT). 

For georeferencing spatial data to a map base, Arc/Info uses an affine equation 
which changes X,Y coordinates and pixel size, thereby losing a linear one-to-one 
correspondence in the process. Arc/Info is a GIS software package and does not 
include a polynomial equation to maintain pixel location and spectral resolution of 
remotely sensed image data. Pixel values can be changed locally. For this reason, the 
use of Arc/Info for georeferencing and image registration is not recommended for 
remote sensing image data. 

Remote sensing data requires remote sensing software for georeferencing. 
ERDAS/Imagine was used for georeferencing. After locating corresponding points 
on an image and on a map (e.g. road intersection), remote sensing professionals 
require software to georectify image data by the following parallel steps: 1) calculate a 
transformation matrix (coefficients) from distance between source ground control 
point (GCP) locations and new GCP location using a least-squares polynomial 
equation (Schowengerdt) to order the transformation (convert coordinates) for a 
linear pixel shift or non-linear pixel bend or warp and 2) resample data which 
assigns image pixels to the resultant grid with options to extrapolate pixel values 
(maintain or change pixel values as required by different types of image processing). 
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Each TM composite scene was georectified to United States Geological Survey base 
maps of coastlines and roads (DLG's) already contained in the State's GIS. A 
resampling method of nearest neighbor was used since the data was ordered in the 
Space Oblique Mercator (SOM) projection with a cubic convolution resampling 
method (and had been convolved once already). The unusual presence of visible 
mountain roads in the summit areas of Mauna Kea and Mauna Loa provided an 
excellent source of control point data to help correct for terrain displacement. 
Including DEM data in the layer stack is another source (Schowengerdt). Map area 
measurements and location control for the raw data mosaic product is based on 
6(32. Om x 32.0m) grid cells and vegetation and other land cover map products are 
based on a single thematic layer of 32.0m x 32.0m grid cells. 

The Hawaiian Archipelago lies in UTM Zone 4, except for the island of Hawaii 
which lies mostly in UTM Zone 5. All geographic databases of the State of Hawaii 
Geographic Information System are in UTM Zone 4. Each of the four scenes in the 
Landsat TM mosaic were georectified to UTM Zone 4, Clarke 1866 Ellipsoid, Old 
Hawaiian Datum, to comply with the State of Hawaii database design which allows 
plotting and printing of the whole archipelago in one UTM zone. The mosaic 
section of this report explains why a model was necessary to process the four scenes 
together into a whole island mosaic. ^ 

There is a need to be aware of the different ellipsoid and datum systems used in 
Hawaii and to know what sources a project is using. Problems occur when sources 
are not known. This affected the project in two areas: change detection image-to- 
image registration with Landsat TM and 1977 Landsat MSS data; and using a version 
of distribution maps of Hawaiian vegetation (hereinafter called "Jacobiveg") Level II 
overlays, to compare Level II nomenclature and ecozone boundaries with the 
NASA projects GPS field site data and ecozones in the imagery. 

The Clarke 1866 ellipsoid is used in the State GIS database and is centered island-by- 
island. The U.S.G.S. uses Clarke 1866 for DLG’s, but not for topographic maps. On 
topographic maps, U.S.G.S. uses the International 1909 ellipsoid and NAD27 Datum 
(Ev Wingert 1995). This was used to georeference the 1977 Landsat MSS data to UTM 
Zone 5 in the previous project. 

The different datums are Old Hawaiian Datum, NAD27, and NAD 83. Old Hawaiian 
Datum is used in the State of Hawaii database design and has four different 
centerings, centered on a latitude / longitude system. Old Hawaiian is sometimes 
mistakenly referred to as NAD 27 which is globe centered (Ev Wingert). 

Problems occur when transferring field data plotted on U.S.G.S. quads into the State 
GIS, if the sources are not known. The ecozone boundaries of Jacobiveg Level II 
digital overlays were displayed offset from the image data in some areas. The source 
of Jacobiveg mylars is stereo pair photos interpreted and plotted onto paper source 
orthophoto quads with the fitting in UTM Zone 5, International 1909 ellipsoid. The 
project source is digital in UTM Zone 4, but a X,Y shift is required on many quads. 
They need to be re projected for the State GIS database to date. 
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Spectral Image Classification 

Digital image classification is the process of deriving land cover and land use classes 
from multispectral images. Through this process, the multiple spectral bands of the 
image are reduced to a single layer of information for interpretation in land 
management or research use. Image classification is a method of processing spectral 
ata rom a satellite or aerial reconnaissance view, which is combined with ground 
view information or field spectra data for mapping and analysis. 

The Landsat TM satellite view is unique, recording radiation bands of intensity 
measurements in visible, infrared, and thermal wavelengths. During the 
classification process, the multi-band image is grouped into a analyst-defined 
number of clusters, based on the intensity values of the image pixels (grid cell). Each 
image pixel is described by a vector through the total number of bands in the image 
file (n-dimensions), or a subset of bands given to the classifier. With spectral data a 
statistical model is used to set the natural variability of vegetation and other land 
cover classes and the classes average. A decision rule determines which class 
vegetation goes to. Statistical decision rules assign the value of one of the clusters to 
the 6-d measurement vector for each pixel in the image. Descriptive statistics are 
saved to a numerical signature file for analysis and different band combinations of 
clusters can be displayed in 2-d graphical plots of histograms or ellipses'^The 
resultant thematic image is colored and interpreted one category at a time, while 
using ancillary information such as photography and vegetation keys, GPS points, 
held sampling information, and vegetation overlays. 

The spectral image analysis process requires knowledge of which spectral regions 
identify different types and level of land cover information needed. The resolution 
omain of the sensor, knowledge of the analysis process and the ground, help to 
locate which spectral regions identify ecozone patterns for the level of land cover 
information needed. 

This means the resolution domain of a sensor and data need to correspond to the 
domain of the ecosystems being mapped. Natural variables of each land cover type 
in a ecosystem need to be identified when producing a classification with 
hierarchical layers, ranging from general to specific. The natural variables of each 
and cover ecozone or transition zone still need to be identified when producing a 
classification with the variables grouped into more general classes. Land cover 
variables improve classification accuracy, because spectral space will go somewhere 
if a land cover is not accounted for. 

Classification Methods 

There are two methods of image classification, supervised and unsupervised. The 
value of a particular method and decision rule depends on how accurately the 
image data is represented for the classification level or change detection that is 
nee e . ^ speed of the program, time for field data collection and interpretation 
are critical factors. Data compression techniques prior to the classification process, 
reduce the dimensionality of information in huge data sets like Landsat TM or 
hyper spectral data, or data from multiple remote sensing platforms 
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An unsupervised clustering is a semi-automated approach to classification, whereby 
the image data points (pixels) are grouped into non-overlapping clusters, based on 
measuring similar pixel values in the image file. The analyst specifies the number 
of clusters, the number of iterations or passes through the data, and a convergence 
value, which is the percentage of data that is successfully assigned to a class before 
clustering (or partitioning) is finished. The classifier assigns the value of one of the 
clusters to each pixel in the image. The resultant clustered image is a gray scale and 
the classes produced are spectral classes. They are interpreted one category at a time, 
while labeling clusters with color and land cover names. 

Although general, land cover derived from an unsupervised classification with 30- 
50 clusters can be processed and interpreted in 2-4 days time, depending on available 
ancillary information, knowledge of the ground, and plotting and printing 
capability. 

In a supervised classification, the analyst manually selects training sites on the 
image monitor for each known land cover in a geographic area. Each training site 
represents a sample of pixel values in the image, which corresponds to a field site of 
land cover on the ground. The values of pixels in the training site are used to 
compute cluster statistics of each band for a corresponding class. The classifier 
assigns the value of training site statistics to each pixel in the image. The training 
site is located by knowledge of the ground with ancillary data such as photography, 
GPS, vegetation and soils overlays, elevation, or other types of field data. The size of 
the sample needed varies with the sensor used and the number of dimensions in 
the image set. The type of classes that are produced are trained information classes 
instead of unsupervised spectral classes. 

In either classification method, the analyst interprets the resultant thematic image, 
with the cluster numbers ranging from low to high, based on the associated image 
pixel brightness values. A cluster number for recent pahoehoe lava is low with 
associated dark pixel values, while cloud cover has high cluster numbers associated 
with high pixel values. In either method, there are tradeoffs of time and method of 
sampling selection. The supervised classification method is more interactive and 
time consuming, followed by unsupervised classification with a large number of 
clusters if a prior knowledge or photography are missing. 

Island of Hawaii - Discussion of Landsat TM Image Classification 

The classification methods described were used for land cover mapping of the Island 
of Hawaii, using Landsat TM data and ERDAS/Imagine v 8.1 which has 
interchangeable classification tools. In this project, four georectified TM scenes were 
classified individually to produce map products. Because the spectral and temporal 
resolutions of the project data are very different, their interaction for land cover 
classification was a concern from the initial Hawaii Island planning meetings in 
April, 1994. Images from two different satellite orbits with different atmospheric 
conditions, one an El Nino and the other a drought year, were going to be part of the 
same land cover classification map product. 
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Additional processing was needed before mosaicing the four scenes to extract the 
temporal land cover information processing for different moisture levels and 
growth states of vegetation, primary succession on lava flows, substrate changes in 
agriculture and volcanically active areas, and atmospheric conditions. 

Spectral Band Combinations Used for Classification and Analysis 

Various band combinations were used to discern the complexity of azonal substrate 
regions and vegetation-substrate cover observed in zones. TM bands 7-5-3 were used 
for analyzing the separation of lava, soil, vegetation moisture, and vegetation 
pigmentation concentrations. These bands also detect the heat from the 'skylights' 
of active lava flows and older, cooled lava. Bands 7-5-4 were used to visualize 
moisture in grasses, biomass, and atmospheric effects along with 7-5-1. Bands 7-4-3 
provided maximum separation (or decorrelation) of soil, lava and vegetation 
biomass. Band combinations 4-5-3 and 4-7-3 were especially useful for biomass and 
maximum separation of Ohia-Koa from other species, i.e. Dense Myrica faya and 
Tropical Ash. However, several vegetation types in the project area displayed dense 
biomass or different moisture contents from the El Nino and drought years, 
appearing saturated in one or the other band combinations. More grass and shrub 
were visible in the vegetation mosaics on the drought dates than on the El Nino 
dates, probably -due to moisture stress and wilting. 

Unsupervised Classification Process - Island of Hawaii 

The unsupervised classification algorithms in ERDAS Imagine v 8.1 were used 
with the four Landsat TM scenes to derive land cover and land use information, 
and to monitor change. ISODATA (Iterative Self-Organizing Data Analysis 
Technique, Tou and Gonzalez 1974) with the minimum distance decision rule was 
used as well as RGB partitioning (ERDAS, 1994). 

ISODATA clustering is a semi-automated classifier, allowing the analyst some 
control by specifying the number of clusters, iterations, and a convergence factor, 
which is the per cent of data classified before clustering is finished. Through this 
process multiband spectral images are grouped into a analyst-specified number of 
clusters, based on spectral groupings of image pixel values around their average. 
The minimum distance statistical rule is used to assign the value of one of the 
clusters to each pixel in the image. Final clusters are represented by a normal 
distribution. 

Instead of using a specified classification system, unsupervised classification relies 
on defining the spectral classes during the interpretation process. The clustered 
image that is produced is a single clustered layer (thematic layer), and interpreted by 
labeling each cluster with color and a vegetation, land cover , or land use name. 

The project team used the ISODATA method so all six TM bands could be processed 
and the number of classes and iterations were varied to distinguish spectral classes 
of land cover and land use. Band 6, a thermal, was left out because of its coarse 120m 
spatial resolution. The algorithm is iterative and cluster means are re-calculated 
each iteration before final clusters are represented by a normal distribution (bell 
shaped curve of probability distribution). 
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Clustering by specifying a large number of clusters helps to represent the full range 
of pixel intensity values in the image data (0 -255). The resultant clusters are 
interpreted, given a land cover name, and grouped together forming spectral 
patterns and land cover classes. Clustering with this approach was used to specify 
200 to 255 categories, and 10 and 100 iterations, using the default convergence to end 
clustering. Analysis requires apriori knowledge of the ground and knowledge of 
how to locate the spectral clusters that identify the land cover classes needed. 
Sample areas from these maps were useful in the field. 

Clustering with 40 classes and 35 iterations per scene was easy to interpret and 
provided adequate separation of classes to provide preliminary map products to 
begin the field survey, although land cover patterns were generalized. The whole 
classification and interpretation process took 1.5 man days to complete each time the 
classifier was processed. These maps were very informative for collecting field data 
to use in the subsequent supervised classification and to monitor spectral areas 
difficult to classify in the supervised classification. 

Two obvious areas needing monitoring occurred where there are two classes that 
are close to each other spectrally located between 5,000-6,000 feet in two different wet 
zones; one area is rainforest with Montane open Koa-grassland infested with dense 
banana-poka, while the other area has Montane Grazed Grassland infested^ith 
Gorse in the mountain pastures. 

The second confusion occurred in the transition zones of mixed grazed grasslands in 
leeward areas where grass-cactus-shrub confused with Montane Wet Eucalyptus in 
Laupahoehoe. During the subsequent field surveys, the area was monitored and 
GPS and non-GPS points were taken to determine their spectral regions. 

Examples of RGB clustering were made during training. RGB clustering is more 
automated and allows for only three spectral bands, which does not take advantage 
of the multi-dimensions of Landsat TM or hyper spectral data. RGB clustering 
performs a principal component analysis prior to partitioning the data one time for 
clustering. The classification is based on systematic sampling of image pixels that 
may pass over spectral regions needed for a land cover classification. 

Ancillary information for the unsupervised process included ground information 
from DOF AW/Hilo, orthophoto quads, the previous MSS classification, a DOFAW 
forest atlas, color infrared photography from Air Survey Hawaii and the ER-2 (U-2) 
from NASA at Ames Research Center. 

Field Work 


Following the generation of the unsupervised classifications for the north scenes, an 
Arc/Info Point coverage was generated to overlay on the unsupervised classification 
map along with roads and management layers. 

A field survey collected data to evaluate the classification of the north west scene 
and to collect ground location data with the GPS for the training sites needed for the 
subsequent supervised classification. While resource managers use aerial 
photography to plot ecozone boundaries and monitor land, they are more familiar 


25 



with a horizontal view of pattern recognition than a reconnaissance view. While 
gathering field data, adjusting to the reconnaissance view requires consideration the 
image date, spectral wavelengths, and the minimum spatial resolution and map 
unit of the land cover-vegetation map product. Relating field ecozone information 
to the spectral regions in the imagery is a critical part of the classification process. 
GPS points provide accurate ground location coordinates of vegetation and other 
land cover in the landscape, for the otherwise difficult job of locating corresponding 
areas on the imagery. 

Initially, field sites were chosen that typified the major wildland vegetation and 
other land cover types to be used for the supervised classification of the two north 
scenes. A total of 41 points were taken over 3 days. A sample raw data image and 
unsupervised classification map were used to avoid taking GPS readings in areas of 
the scene containing clouds and shadows, but four points were taken near these 
areas because they provided needed information for the map. These points were 
saved for future use, but could not be used with this data. For each point, a site 
description and photographs were taken. GPS location information for the major 
land cover and land use was used in the supervised classification. 

Additional GPS points were collected in field sites, while assessing the accuracy of 
the supervised classifications of the north scenes. The goal was to locate land Tover 
areas needing correction, while collecting new field samples. Points were taken to 
locate transition boundaries of some ecozones seen on the imagery and misclassified 
on the maps of the north scenes. 

The abundant transition zones were found to contain spectral regions representing 
diverse land cover type variables. Examples are two dominant grassland areas of 
Parker Ranch, mixing with two other grass types and shrub types that change based 
on elevation, moisture-climate, and grazing management patterns. (Another 
example is the re vegetated lava class described in the Island of Hawaii-Classification 
System section). The same occurs for all natural variability found in vegetation 
ecozones and patterns at the landscape level. 

GPS Methodolog y 

The field team used a Trimble 6 channel Pathfinder Basic with an external antenna 
for obtaining GPS position data. GPS uses the L-band group of radio frequencies to 
collect signals locating ground position, and documents these positions in the field. 
(Data is collected using WGB84 frequencies). Field data were differentially corrected 
to the local base station on the big island. The process compares your field position 
data against stationary receiver data collected by a base station. The base station at the 
National Historic Park, (National Park Service) of Kaloko-Honokohau in the 
Kailua-Kona area of the big island was used to improve accuracy. 

The corrected GPS points were projected into UTM Zone 4 to include the whole 
island with the State of Hawaii database design. Accuracy of approximately lm is 
available with differentially corrected GPS points. Differential correction of the GPS 
data was done by DLNR staff using PFINDER SOFTWARE version 2.4. GPS data 
processing yielded text files with UTM coordinates. The text files were used to 
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generate a point coverage in Arc/Info., and the collected field information was 
entered as attribute data for each point. 

Non-GPS Points 

Field samples of areas lacking GPS points were plotted on maps and raw data prints 
throughout the project. Hamakua Sugar Co. and Kau C and H Sugar Co., provided 
small and large scale maps containing field numbers for sugar and macadamia 
crops, and field data of the dates of the sugar crop planting cycle, revealing the end of 
the sugar industry on the big island. 

Raw data prints and classification maps from an 8.5Xllin. printer, photography, and 
topographic maps were used to plot field samples of some of the differences found 
in the main ecozone domains of Hawaii Volcano National Park, Puna, Kona, and 
Kahala areas. 

APPENDIX B - Ancillary Field Data 
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Supervised Classification - Island of Hawaii 

The supervised classification programs in ERDAS/Imagine v 8.1 were used to 
process the four multiband TM scenes to derive a single, corresponding layer of land 
cover categories for the Island of Hawaii. The supervised classification for the Island 
of Hawaii began with the north scenes which contain the Hamakua coast, and is the 
initial project site for processing Landsat TM and the previous projects MSS data set. 
In the supervised classification process, a numerical file of spectral signatures was 
built for use with the maximum likelihood decision rule. The supervised 
classification process involves three areas: training the classifier, classification with a 
decision rule, and product evaluation. 

Training is the process of selecting samples of image data that represent ground land 
cover to develop numerical training signatures that guide the classification of the 
whole image. The training sites are used to build the numerical signature file of 
spectral signatures for processing each land cover information class needed on the 
map. Regions on the image represent spectral data samples of the field sites visited 
on the ground, or found on photography. The training process is usually repeated 
before the desired signatures are produced. Descriptive statistics are displayed for 
training sites of each land cover type and for combinations of land cover types, that 
may confuse. During the training process, areas are flagged for additional field Work. 
The unsupervised classifications are informative visual tools during this process. 

The classifier computes spectral patterns in the image data training set using a pre- 
defined statistical decision rule. Each pixel in the image is compared numerically to 
each training site category in the numerical signature file and categorized to the land 
cover class it is most like. The category label assigned to each pixel in this process is 
recorded in a corresponding cell of the classified image. 

The supervised classification image is evaluated on a monitor, tabular, or map 
forms. Following evaluation, training sites are added, modified, or deleted to derive 
the data representation of each field site that corresponds to a land cover or land use 
site found on the ground. 

Processing for Subtropical Differences - Island of Hawaii 

The Hawaiian Archipelago has physical, biological, and manmade ecozone 
processing requirements that are different from continental environments. 
Contributing physical ecosystems are vast amounts of water and cloud cover,^coral 
and anchialine pools along the coast, small land areas with extreme elevation 
changes, on-going volcanic eruptions with skylights of molten lava from networks 
of lava tubes, and both substrate and atmosphere contain VOG. Contributing 
biological ecosystems change from tropical cloud forests to woodlands, sub alpine 
shrubland, and alpine cinder fields and deserts. Contributing manmade ecosystem 
include quarries and gravel pits of lava, fish ponds, tropical agriculture of sugar 
cane, abandoned sugar fields, macadamia, coffee, and nurseries with tropical trees 
and flowers. The need for multiple scenes to cover a relatively small area is also 
common. 
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The multitemporal data sets for this project, increase the dimensionality of the 
imagery to process, by requiring the addition of spectral sub-classes for different 
dates of available field sites. In the future, each new set of image data will have 
unique characteristics that also require some custom processing. 

Developing a spectral signature file for the island of Hawaii began with appending 
water and cloud signatures from the files generated by the four unsupervised 
classifications to account for orbital and atmospheric differences of the adjacent 
scenes. A signature file for the north west scene was used as the base for all four 
scenes of the big island since this scene contains most of the island land mass and 
the test site. The water and cloud signatures were used during the classification 
process and the water signatures were also used in a masked classify application. 
These techniques are especially helpful for image processing in Hawaii's complex 
environment. 

Water masks were made in ERDAS/Imagine and Arclnfo to exclude land cover 
information from processing and to combine layers processed separately. Water 
masks were made at the coastline and to include off shore islands. The masks can be 
used in subsequent classifications of this image or registered to new data. In 
Imagine, the combined water classes from the four unsupervised classifications 
were used as a mask, as were the infrared bands to determine the water-Wnd 
boundary. The water classes from the unsupervised classification were easier to use 
than the infrared bands of the different dates. 

For clouds, determining limits for the brightness values of cloud cover, urban, and 
dry grass was also required to avoid misclassification of cloud edges as dry grass. (A 
post-classification mask was made for the Cloud to Urban category, which is the last 
category in the masked classification mosaic). 

Processing for Subtropical Manmade Ecosystems-Land Use 

The supervised classification training process started with the sugar plantation data 
from Hamakua Sugar Company in Paauilo, including the Kahala area which ended 
sugar production in 1977, the same year of Landsat MSS data. Ancillary data 
included sugar maps with field numbers, and data for crop calendar and the 
termination of production. Color infrared photography was also used for 
interpretation and to monitor land use change. Individual training sites were 
developed with the tools in Imagine and edited in the signature file. 

Abandoned cane and pastures with Christmas berry, Ironwood, Koa haole, or kiawe 
were found on the ground in north Kohala and the Waipio side of Hamakua, but 
the imagery was covered with clouds. The increase of nursery rooftops was noted in 
Hilo and Puna areas, along with macadamia and papaya fields on lava or grass. 

Land use plantations of sugar and mac nut were resolved in Lowland areas 
containing large fields. The small diversified ag fields containing mac nut and other 
crops on grass need multiple resolutions of image data or photography, field work, 
and GPS to be resolved. The small Lowland Diversified Ag fields of ca. 1-5 acres and 
Montane zone mac nut orchards confused with some tropical wildland vegetation 
and gorse, and were eliminated from the signature file. 
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Some land use conversion was observed during field survey of Waikea and Puna 
areas. Small nurseries containing palm confused with wildland palm and were 
grouped into the Lowland Mixed Forest class. Papaya was not trained on because of 
the need for historical crop calendar information. Urban areas and nursery rooftops 
confused with cloud and were masked in the densest areas as the last information 
class on the masked classification. 

Processing for Subtropical Biological and Physical Ecosystems 

The GPS points were used to develop signatures by locating image data points with 
corresponding ground field sites of wildland vegetation and physical ecosystems. 
The GPS points and their attributes were displayed over the TM imagery to locate 
the corresponding field site boundaries. GPS and non-GPS points were used with 
other field data information to locate a cluster of points representing the land cover 
ecozone where they taken. Spectral and spatial response patterns of the different 
vegetation and land cover were analyzed for natural variation. On the image data, 
the GPS point correspond to a vector through 6-dimensions of the image coordinate 
(for TM bands 1-2-3-4-5 -7). The GPS points were especially useful in the rain forest 
areas with dense biomass; in some locations of closed canopy GPS positioning 
required moving to a slightly open area to collect GPS points under forest canopy. 

The north scenes contain five forest classes with banana poka in Hamakua and 
various band combinations reveal pigmentation and moisture stress levels. The 
project’s GPS points were used with Jacobiveg to compare the new GPS point data 
with Jacobiveg Level II polygons. Coarse distribution maps were also used for 
support information (F. Warshauer, J. Jacobi, A. La Rosa, J. Scott, Cliff Smith, 1983). 

Jacobiveg was used to see how the training signatures with GPS points were 
classifying at the canopy level on the Kalapana and South Point scenes that lacked 
ground truth initially. The vegetation ecozone coverage of Jacobiveg Level II were 
also used with the GPS points and the classification summary of W. Gagne and L. 
Cuddihy Pratt (in Vegetation, Manual of Flowering Plants, 1994). Differences in 
ecozone boundaries of raw imagery and the classification map were noted. 

Training sites were generally located on the raw data using available ancillary 
information including color prints from Imagine viewer annotated with field data, 
GPS coordinates, aerial photographs, field data supported by topographic maps, 
orthophoto quads, an outdated forest type atlas, sugar plantation maps,. and 
vegetation overlays. 

During the analysis process, several processing techniques were developed to 
visualize some of the impact caused by the temporal, seasonal, and elevation 
differences. The signature file developed for the north west scene (El Nino date), 
was used to process the north east, Kalapana, and South Point scenes. Initially, 
multi-date signatures were not developed to visualize the effects and a mosaic was 
built. In Imagine, the mosaic output is a gray scale image associated with a gray scale 
bar, until the image is loaded as pseudo color with the color bar or the R-G-B values 
from the original classified images. 
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Several land cover types did not correlate across scene boundaries due to temporal 
and seasonal differences of some land cover types. This caused the same land cover 
and geographic area to classify in different information class positions of the 
classified image and gray scale bar. The color look-up-tables (LUT) from each scene 
were loaded manually in the mosaic and flooded the adjacent scenes, where 
classified as the same land cover type but a different information class number or 
where classified as a similar land cover type. Part of the problem is a raster hardware 
and software constraint, lacking the ability to trap class color boundaries the way a 
GIS package is able to. 

Once the temporal differences were resolved, several land cover types still did not 
correlate across scene boundaries due to a few signatures being unique to one or two 
scenes, so the number of information classes were different for corresponding 
scenes. This also caused the same land cover and geographic area to classify in 
different information class positions of the classified gray scale image and gray scale 
bar. 

On the north west scene, the spectral confusion between grass-shrub-cactus and 
eucalyptus was resolved by developing spectral signatures with vegetation change 
based on elevation and moisture differences. Five information classes were 
discerned for the mixed grass-shrub types in the transition zone. These five classes 
are colored with Fountain grass-Mixed grasses on the map, but they can be colored 
individually in the raster or vector files. Classes range from Buffelgrass-Fountain 
grass, to Buffelgrass-Fountain grass-Guinney-Kukuyu grass with shrubs. 

Multi-Date Classification for Whole Island Mosaic 

The final solution to mosaic the whole island map, was to construct one main 
signature file, which was used to classify each scene individually. The main file was 
modified for each scene by weighting the land cover signatures unique to a 
geographic area of the island or to one scene. This method requires an inordinately 
large number of repetitive training sites to develop a higher dimension of statistics 
for the temporal dates, and the information is repeated in the class names list of the 
raster attribute table of the final classified mosaic. 

Training site samples were chosen that uniquely represented the date of their 
corresponding land cover and land use field site. This means a higher dimension of 
statistics was needed to account for the visible moisture and biomass levels_, and 
growth states. The multitemporal data represents a vast amount of added 
information and processing complexity needed to extract information. Training sites 
were evaluated prior to classification by image and numerical analysis. Histogram 
displays of the descriptive statistics in training sites, were used to insure separability 
and determine how the data covary. Descriptive statistics of pixel values for each 
training site were used to identify the spectral signature that best represented a 
corresponding land cover/land use class. Multivariate analyses of cluster plots and 
feature space images were used to quantify the separability or spectral difference 
between training sites and to anticipate misclassification of classes. 
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Pixels in each image are assigned to a class by the multi-date maximum likelihood 
decision rule unless a mask is applied to omit a specific value or range of values, 
representing land cover. 

On the mosaic, there are vegetation types with only the core area of their 
concentration mapped to a intermediate class. More field data and processing is 
needed to map their extent with spectral regions for land cover class variables. 
Examples are Ash, Sugi, Eucalyptus, and Gorse which appeared on earlier maps of 
the two north scenes, and Myrica faya-Ohia. 

Dense Myrica faya (fayatree or firetree) with a closed canopy was mapped in the 
Hawaii Volcano National Park area and in Hamakua using coarse distribution maps 
to locate field sites (Tunison and Camrath 1992; and Whiteaker and Gardner 1985). 
Fayatree occurs in two mixed classes of the classnames list/Legend in Appendix C. 

Additional data was collected but not processed, and more data is still needed to 
date. Scattered pixels of Gorse on the map were field checked and do occur north of 
the Saddle Rd. It is unknown to date what signatures could be developed for 
intermediate classes and how these would classify, but additional classes exist. 
Renaming existing classes to Myrica faya on the individual scenes is possible, but 
changing the color with the current mosaic constraints would cause the class to 
appear over extended on adjacent scenes. This area could be masked on the mosaic, 
but the color tables for the corresponding Ohia class number on adjacent scenes 
without Myrica faya would be corrupted or change to the color of Myrica faya or 
Myrica faya-Ohia. 

Signature development for these land cover types was successful on the individual 
scenes, and in decision rules of the intermediate iterations in the main signature 
file. The signature was minimized in the last iterations of all four scenes when new 
field data was added. This is fixable through successive image processing. The cell 
array information in Imagine documents the current iteration. 

Visual systems are needed on all islands where land managers and researchers are 
left without visual tools for the lands they are responsible for. A portable system 
would include outer island resource personnel in the processing, and help everyone 
involved. 

APPENDIX C - Classification System 


Whole Island Mosaics 

The mosaic process assembles the four adjacent scenes into the whole island 
together with Maui. This process requires georectification of the scenes to the same 
map projection (UTM, UTM Zone, spheroid, and datum) prior to the mosaic process 
for consistency with State of Hawaii database design. The process was repeated for 
the four georectified composite scenes and the 4 supervised classification images. 
Initially, all of the automated menu options were tried but these placed conspicuous 
and unacceptable lines or wedges through Mauna Kea, the unfortunate location 
decided by the automated menus. The adjacent scenes were of the extreme temporal 
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differences of a El Nino year (June 1, 1990) and a drought year (October 11, 1991), and 
from different seasons; differences in the moisture level and chlorophyll 
concentrations of adjacent scenes were quite obvious; for example, biomass of 
understudy and ground cover was more apparent in 'local area' mosaics of open 
vegetation. 

Clouds in the 11% overlap area between scenes were cut off leaving a sharp 
boundary line where there wasn't cloud cover on the adjacent scene. Custom 
models were built in order to have flexibility over where to place the mosaic 
boundary and to reduce cloud cover in the overlap area, thereby recovering cloud 
free data in the adjacent scene. This process recovered image data of some of the 
forest reserves in the overlap area between adjacent scenes. 

The process works only in the 11% overlap area of adjacent scenes because the area 
of interest (AOI) processing used in the modeling requires the geographic area of 
second image to fill in data for the current blank area. If there is cloud on the first 
scene, the grouped AOI's point the model to the same geographic area in the second 
scene of the model, which fills the overlap area of the adjacent scene with data. 

The final step in the image classification was to mosaic the classified images into a 
composite consisting of the whole island of Hawaii. This process was Jime 
consuming because of the temporal differences in the data (adjacent El Nino and 
drought years). The classification and mosaic process of each of the four scenes was 
repeated when additional ground truth was collected. 

Importing the classification RGB color tables, class names column, and the addition 
of a legend column were required to group the individual class names in the mosaic 
image (Appendix C). A parameter to automatically build this information from each 
of the classified images and a trap for boundaries is needed in the raster processing. 

APPENDIX D - Mosaic Process 


MAP COMPOSITION AND MASKING 

Masking is a type of area of interest processing used to select land cover information 
for processing, combine layers processed separately, or edit misclassification. This is 
especially helpful for processing in Hawaii's complex environment. The various 
masking tools in Imagine and Arc /Info were used. _ 

Masking clouds was used in the initial part of the whole island mosaic process. 
Masking was required to assemble the four scenes into an image containing the 
whole island, and to exemplify renaming of misclassified land cover types on the 
supervised classification land cover map. See the section titled Mosaic. 

The whole island mosaics were processed in Imagine v8.1 on the HP platform, 
which required building a graphical model since the menus did not work 
completely. Areas were grouped in multiple or nested polygons, extracting only the 
areas needed to assemble the four scenes together. 
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Water masks were difficult to make using the infrared bands of the multitemporal 
scenes and mosaic, so signatures for water classes were extracted from the 
unsupervised classifications and appended to main signature file instead. They can 
be used this way or to classify water and land separately to improve per class 
accuracy. 

Examples of masking the classification using the Spatial Modler were processed on a 
temporary license of ERDAS/Imagine v 8.2 for the SUN Microsystems platform at 
The University of Hawaii Geography Department. Urban was intentionally added at 
the end of the image processing to avoid confusion with other land cover types in 
the classification. The software menus for recoding classes with new class numbers 
was used with a area-of-interest polygons digitized on the color monitor. This 
streamlined processing considerably and was somewhat of an inter-active process. 

The raw data color composite of bands 7-5-3 and the masked land cover map were 
exported from Imagine and into the Arclnfo graphics files for map composition. The 
map production AML’s contained in the State of Hawaii /GIS were used to compose 
these maps. 

APPENDIX G - Map Composition 


Accuracy 

Accuracy assessment consisted of several field surveys for map validation, 
renaming, and field data collection. The land cover map produced at 1:250,000 and 
1:24,000 scale is an overview level or base where field sampling of ecozones is 
unprocessed or missing. This reflects the domain of the sensors resolutions 
combined with the lack of funds for the level of land cover variables present in 
Hawaii ecosystems and temporal data (i.e.. cloud cover). 

If the image classes are renamed to a coarser level where needed, classification 
accuracy improves. If the imagery is classified on a 7.5’ quad basis, and combined 
with GPS, or multiple sensor platform levels (and georeferenced scales), and 
detailed vegetation or soils overlays, the level of spectral classification accuracy 
increases. 

Additionally, some recoding of classes or masking can be done on the GIS side (with 
reclass) to improve classification accuracy before dissolving polygon boundaries, into 
larger map units. Both are typical parts of the project flow. (See Vectorization). 


Change Detection Image 

Image data that is georeferenced together from multiple sensors, provide multiple 
resolutions and additional dimensions of land cover. If the data are processed and 
plotted at a specific scale, then multiple scales of land cover patterns are available. 

The change detection image provides greater spectral separability between the 1977 
Landsat MSS mosaic and the new multi-temporal Landsat TM mosaic. Information 
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is visible in this image that may not be obvious in the MSS or TM by themselves. 
Registering new data from different sensors and resolutions to the existing image 
database allows processing and plotting at multiple scales to produce multiple scales 
of land cover patterns (i.e.. under flight). Processing these data sets requires 
historical information in conjunction with current field sampling. 

The 1977 MSS data from the previous NASA project was a mosaic assembled of two 
scenes from one satellite path (N-S orbit) and from two rows, with analog pixels of 
79 X 79m digitized to 56m X 79m pixels of the ground. The data were quantified to 6 
bits ranging from 0-63. The scene dimensions were 2340 lines X 3240 columns 
covering 185 X178km on the ground. The mosaic was processed during the previous 
NASA project (C. Hogan, G. Fosnight, Ed Petteys, C. Tasaka, D. Morse, P. Costales, M. 
Buck, and Pat Chavez; principal investigator Ken Nishioka). 

Pat Chavez of USGS, Flagstaff, Arizona generated spectral transformations on the 
1977 MSS mosaic during the previous NASA project, since he had developed 
algorithms which NASA/ Ames did not have at the time. (The mosaic from the 
previous NASA project appears in Volcanoes of the National Parks - MacDonald, 
Wright, and Erickson; 1989). 

During this project, the MSS data was transferred from 9-track tape to CD in UNIX 
in BSQ format, using the UNIX command 'tar' at University of Hawaii PG/HIG. 
The MSS data was transferred again from CD to the file server at State of Hawaii GIS 
lab by systems analysts in secured environment. 

The ERDAS/ Imagine Import utility was used to load each band into Imagine with 
generic binary read option, to obtain the ERDAS image file format (.img file format). 
A graphical stack model was used to make a four band color composite of the 1977 
MSS mosaic. 

An image-to-image registration of 1977 MSS big island mosaic and the new TM 
mosaic, were combined producing a new ten band image. In the previous project, 
the 1977 MSS mosaic was georeferenced to UTM Zone 5, using topographic maps. 
New control points were located on the images and the 1977 MSS was registered to 
the Landsat TM in UTM Zone 4, Clarke 1866, Old Hawaiian Datum. The Clarke 1866 
ellipsoid is centered island by island. On topographic maps, U.S.G.S. uses the 
International 1909 ellipsoid. 

False color or- 'pseudo-color' composites of different band combinations -were 
displayed in which the reflected ground features do not always appear in natural 
colors. The order of the bands in the composite correspond to the red, green, and 
blue (RGB) color guns of the monitor. Band combinations were charted and areas of 
change in vegetation composition, biomass, moisture, atmospheric effects, and new 
lava flows were observed, often times looking opaque. 

Jacobiveg Level II vegetation overlays for five quads were displayed over areas of 
change with GPS field points or non-GPS field points. The most obvious areas of 
change are cloud cover, reduction in the number of sugar fields, moisture 
conditions over the Hamakua and Volcano area forests and grasslands on the 
drought dates, Hilo, and the dense biomass increase of the non-native Gorse, in the 
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high pastures East of Mauna Kea. Jacobiveg and the 1977 Landsat MSS data provided 
needed historical information to monitor gorse. 

The ten band image was compressed and deposited on-line at The University of 
Hawaii Geography Department. In the final stages of this project, change detection 
analysis was performed on the combined TM/MSS image file by Renee Louis and 
Eric Yamashita during a Remote Sensing Research Seminar conducted by Dr. 
Matthew McGranaghan. Seven (7) map compositions were completed using a 60-day 
ERDAS IMAGINE version 8.2 license and plotted on an HP650C. 

The following six (6) map compositions were chosen to demonstrate the power and 
versatility of this technology. The map titled Gorse and Koa-Banana reveals the 
extent of pristine forest depredation due to the invasive nature of two exotic plants, 
Gorse and Banana Poka. The map titled Mauna Kea Area Reserve provides 
evidence that Landsat TM bands 7 and 5 are excellent for analyzing the separation of 
lava, soil, and vegetation. The Keyhole map demonstrates the land cover change 
detection capabilities of this technology by showing the keyhole pattern in three (3) 
images, each with different band combinations (MSS data only, TM data only, and 
TM/MSS combined data). In the map titled Lava Flows, the lava flow which began 
in 1983, and is still active today, is pointed out in the image showing both TM and 
MSS data. The Northwest Hawai'i Golf Courses map clearly depicts the addition a 
new golf course and the modification of existing manmade features. The last map, 
titled Hamakua, Hawai'i Island shows both the combined data and the TM extracted 
classification of this area. 

APPENDIX F - Change Detection 


Vectorization/GIS Layer 

The classification mosaic with areas of example masking was vectorized on a 7.5 
minute quad basis (1:24,000 scale) at the GDSI office in Manoa. This product contains 
1.5gb of data which can be loaded to the file server one quad at a time as needed. The 
Arc /Info command Polygrid was used to generate a 1,000 m buffer into the ocean 
surrounding the big island to include 'the new land' from the eruption of Kilauea's 
Puu O' o and to include off-shore islands. The buffer was adapted from the one used 
on the two island mosaic products, which blank out the water at the coastline. _ 

Classes were renamed or renumbered to a limited extent with recode and mosaic in 
the Imagine v 8.2 because of the LUT's processing constraint of a temporal image 
mosaic. The same problem may exist using the reclass command in Arc, which 
essentially does the same thing prior to dissolving boundaries. Trapping class-color 
boundaries and logic can be used with the reclass command and knowledge of the 
ground, the analysis process, and map composition in a team approach improve the 
process. Care must be taken not to re-introduce areas of previous confusion in the 
classifier. The reclass command does not reclassify the data, since there is no 
quantitative spectral information once the data is in a GIS package. The prime color 
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palette used in Imagine was R-G-B, however Approximate True Colors was also 
used. It is important to use one color palette if the image is going to be vectorized. 
There are places where the color selector wheel truncates the R-G-B values and the 
longer R-G-B values in the cell array need to be exported. This is a problem where 
every pixel is used in the dissolve process. The class names and legend lists are not 
carried over to the mosaic and were imported from the individual classified images. 

Technology Transfer Conclusion 

One of the main features of the project is technology transfer to the State of Hawaii. 
Time was spent developing a preliminary remote sensing classification system to 
use in Hawaii's tropical island environment. Landsat TM proved to be powerful for 
this complex environment and incorporating multiple sensor platforms in key 
areas would be beneficial. The consultants and university taught project participants 
about remote sensing, image processing, and GIS, and the State of Hawaii taught 
about Hawaiian vegetation. Project participants consider the project a success. The 
main project objectives of technology transfer for the production of a digital base 
land cover classification map, raw data island mosaic, and change detection image of 
the TM and 1977 MSS mosaics are met. Project participants put together a multi 
media demonstration using GPS points and field photos with the raw mosaic ©£ the 
big island. 

The products from this project are being used in several ways. Office of State 
Planning is using the LSVIEW program to demonstrate the project to visitors and to 
monitor change in Hamakua land use. Division of Forestry and Wildlife is using 
the map product information in Arc View and will plot the vectorized 7.5 minute 
quads for big island forestry and wildlife as needed. Hilo forestry references the 
various 8.5X11 hard copy prints and map at 1:250,000 since they do not have access to 
a image processing/GIS system. University of Hawaii Geography Department 
graduate students charted their own band combinations and made 8.5 X 11 in. digital 
posters describing land cover changes in the various band combinations and put 
them on the Internet. Other government agencies and commercial companies 
throughout Hawaii are interested in the technology and the products. 

APPENDIX H - Project Products 
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Recommendations 


Building a remote sensing/GIS work environment includes data from different 
sensors to obtain multiple resolutions, aerial photography, stereo plotter for 
viewing and delineating ecozone boundaries, field spectrometer, image data, map 
and data storage cases, hardware/software for image processing, GPS, and GIS, and 
peripheral devices for input and output. 

Input devices include CD reader, various tape media drives, table digitizer, and 
scanners. Output devices include black-and-white and color printers, CD writer, 
plotters, camera and image setter output. 

Cover types with only the core area of concentration mapped need funds for 
personnel to collect and process data. A portable system to take image data and 
processing to resource managers and researchers outer island is needed, along with 
color prints of the data and other information needed to develop more spectral 
classes. 

Register new TM data, and data from different sensors and resolutions to the 
existing image database. Process and produce products at multiple scales to build a 
remote sensing database with multiple scales of land cover patterns for the 
requirements of different land managers and administrators. Incorporate under 
flights, elevation, slope and aspect data to use as a mask if needed, add new multiple 
resolution data sets including radar if necessary. Maintain Jacobiveg vegetation 
overlays by combining GPS point data, image data, and Arc. Use as a mask to classify 
on 7.5 minute quad basis in addition to the base classification provided in a 
classification mosaic. 
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REPORT ON GEOGRAPHIC INFORMATION SYSTEM (GIS) SUPPORT 


Geographic Decision Systems International (GDSI) provided GIS support to the project. This 
support had three major components: 

1) GIS Training 

2) Application Development 

3) On-going GIS Support and Trouble-shooting 

1) GIS TRAINING 

The State of Hawaii has largely standardized on Arc/Info software for their GIS work. Both 
OSP and DLNR have Arc/Info running on the same Hewlett-Packard UNIX workstations on 
which the Imagine software is running. GDSI provided a one week course titled “Introduction to 
Arc/Info” which covered the basics of using Arc/Info for digitizing, editing and displaying 
geographic data. The course included a textbook titled “Understanding GIS, The Arc/Info 
Method” and a complete set of course notes including descriptions of the State GIS and GIS- 
standards. Each student was provided the textbook and course notes to keep after the course was 
complete. Because most of the project participants already had GIS experience, only two 
students took the course, one from OSP and one from DLNR. Since completing the course, both 
have been able to use Arc/Info independently and continue to expand their knowledge of 
Arc/Info and GIS. 

Recommendation: It is important that personnel working with remotely sensed data be aware of 
GIS technology and how to utilize GIS data for image rectification and classification, and how to 
use GIS for display and map preparation Although Imagine has map making tools, the remote 
sensing consultant reported problems with the HP version of Imagine in this regard. Therefore, 
Arc/Info was heavily utilized for map making. In the future, GIS training is necessary not only 
for the technical skills for using Arc/Info, but also for the introduction it gives to spatial data and 
map projections. However, if the map making tools can be used in Imagine, it would probably 
be more efficient to use those tools for map-making rather than using Arc/Info. 

2) APPLICATION DEVELOPMENT 

Part of the remote sensing exercise was the collection of field data for classification. In addition 
to noting the land cover at a site, project participants also used Global Positioning Systems 
(GPS) to locate their position in UTM coordinates that could be later used in Arc/Info and 
Imagine. They also took photographs at many of these sites. Although it was relatively easy to 
prepare paper maps showing the location of the GPS sites, project participants wanted the ability 
to see the photos taken at site to help them remember the site and to assist in classification. 

GDSI used the Arc/Info AML programming language to develop an application to do this. We 
also had to scan in the photographs taken at the field sites. The photos were scanned using an 
8.5" x 14" 24bit color desktop scanner at both 150dpi and 300dpi. It was found that 150dpi was 



sufficient to capture the photo on the screen and was therefore used to save space since these 
files can be rather large. The scanned photos were processed on a PC computer at GDSI, 
sometimes with adjustments for contrast and brightness, and saved as TIFF files. These TIFF 
files were transferred to the UNIX workstations at GDSI across the network. The TIFF files 
were then available for display from Arc/Info. But we needed a way to link the photos with the 
GPS sites and to display the LANDSAT image underneath along with pertinent G1S layers such 
as major and minor roads. Thus the need for the application which we call LSVIEVV. 

LSVIEW is totally menu driven using popup windows on the UNIX workstations It is easy to 
use without requiring extensive experience with either GIS or remote sensing. After starting 
LSVIEW, the user sees the Big Island outline and optionally the LANDSAT mosaic. The user 
can turn on and off other GIS layers by clicking buttons on the menu. The interesting GIS layer 
to display is the GPS layer. When the GPS layer is turned on, small red circles appear on the 
map at all of the GPS sites. The user can then click on a site and, if there is a photo of that site, it 
will be displayed on the screen in a popup window. If there is more than one photo, multiple 
popup windows will appear until all the photos for that site are displayed. Closing the popup 
window removes the photo from the display. 

LSVIEW also allows the user to zoom in on areas of interest for a closeup view. Each area .of 
interest can be named so that it can be immediately returned to later in the session or in 
subsequent sessions. Users can also zoom in on the photos themselves to try and see more 
detail, although for this to be more effective the photos would have to be scanned at a higher 
resolution, perhaps 300 dpi or higher. 

LSVIEW is being used by project personnel and has been demonstrated live on a UNIX 
workstation at the Hawaii State GIS Conference. People are very interested in the mosaic image 
and in how the photos of a site can be popped up on the screen. A copy has also been sent to 
NASA for an earlier demonstration. 

Appendix E contains annotated screen captures to give a quick tour of LSVIEW capabilities. 

Recommendation: LSVIEW made the LANDSAT data more accesible to non-experts and was 
usefijl in providing linked photos to help with classification. There is now software available 
(ArcView2) that could be programmed to have similar capabilities to LSVIEW but with the 
additional advantage that the programs can be am on UNIX, PC and MAC environments. In the 
future efforts should be made to implement LSVIEW capabilities in ArcView rather try to 
enhance the functionality of LSVIEW. Now that a classified image exists, it would also be nice 
to be able to display either the raw mosiac or the classified image. Also, when displaying the 
raw mosaic, it would be nice to choose which bands would be used. 

3) ON-GOING GIS SUPPORT AND TROUBLE-SHOOTING 

Due to the extended time-period for this project, on-going GIS support and trouble-shooting 
became another major paa of GDSI’ s contribution to the project. We stopped in at the lab 
periodically to see how things were going and also worked with OSP project personnel to 
modify their AML programs for map production to generalize them and make them more 



automatic. GDSI was also on call to help whenever there was a problem with Arc/Info and the 
GIS. Although there were many minor troubles, two major areas bear individual attention: 

A) Transfer of Imagine color pallettes into Arc/Info 

Because the plotting part of Imagine didn’t seem to be working on the HP 
platform, it was decided to do the plotting of the maps using Arc/Info. However, 
project personnel were spending lots of time trying to get the colors in Arc/Info to 
match the on-screen colors in Imagine. GDSI worked with OSP personnel to 
develop an AML to automatically transfer the color pallette from Imagine to 
Arc/Info. The user need only save the color pallette as an ASCII file from 
Imagine, then start Arc/Info and run our AML. The AML reads the ASCII file 
and creates the color table and legend key to be used for the Arc/Info map. It all 
works automatically. This AML greatly increased the project personnel’s 
capability to produce maps. 

B) Image conversion for plotting 

Once the map was created in Arc/Info it was necessary to plot it out on a large 
format color plotter. Both GDSI and OSP use the HP650C color inkjet plotters 
which can plot up to E-size in color at 300dpi. However, the map graphics files 
(Encapsulated Postscript) were so large, due to the large size of the raw and 
classified images contained therein, that the plotters ran out of memory before 
completeing the map. The HP650C uses the RTL graphics format internally and 
apparently had to store the entire Postscript file in memory and also compose the 
page to plot in RTL in memory at the same time. 

GDSI had an upgraded version of Arc/Info that the State did not yet have. This 
version had an RTL command that purported to create an RTL format file that 
could be sent direct to the plotter. So we took the graphics file and used the RTL 
command. We still ran out of memory but were able to reduce the size of the 
image until it plotted out. These were the first plots of the raw image that was 
then taken into the field. In the meantime, we continued to investigate why there 
was still an out of memory error. It turns out that the Arc/Info RTL command 
doesn’t produce pure RTL files but rather a combination of HPGL-2 and RTL. 
Therefore, the plotter still had to store the the input data in memory while it 
composed the page in memory. Thus the continued out of memory error. 

An enquiry over the Internet (comp infosystems.gis) brought a majority response 
that a software package calle Image Alchemy would be what we needed. We 
ordered the software ($1000 at the time, now $495 on UNIX, $395 on PC) and 
gave it a try. Image Alchemy composes the RTL graphics page using UNIX 
resources rather trying to do it on the plotter. It can take large amounts of hard 
disk space (700 Mb) and large amounts of time (4-6 hours), but it creates a pure 
RTL file that can be spooled to the plotter and plotted as it is recieved This 
removes the plotter memory limitation. 



All subsequent large format color maps for the project have been produced using 
this method, both at GDS1 and at OSP. An upgrade for Image Alchemy (1.8) has 
come out since we started (with 1.7) and we have noticed that the colors produced 
by 1.7 are lighter and more pleasing than those produced by 1.8. There are 
probably switch settings to compensate for this but for now we used the Alchemy 
1.7 for the final project maps plotted at GDSI. 

Recommendation: The color pallette problem presumably would not have occurred if the maps 
had been produced directly by Imagine, however since our AML has solved the problem, and 
since project personnel are more proficient at producing maps with Arc/Info than with Imagine, 
they will probably not switch back to Imagine even if the map production problem on the HP is 
fixed. This project was done using Imagine version 8. 1 and there is now version 8.2 in which 
the problem may have been fixed. Also, unless Imagine has a pure RTL driver similar to Image 
Alchemy’s, the same problem with memory will occur when plotting. Image Alchemy or similar 
products will continue to be necessary so the ability to solve this problem has been a valuable 
contribution of this project to both remote sensing and GIS users in the State Several other 
Federal, State and County agencies have adopted it after seeing our work on this project. 

CONCLUSION 

There is great interest in the GIS community in Hawaii for using remotely sensed imagery for 
landcover classification and for landcover change analysis. Hawaii’s unique natural heritage and 
wide range of diverse ecosystems makes this a challenge, but one that governmental and private 
organizations are eager to overcome. In this project we have shown that Arc/Info and Imagine 
can co-exist and can work together to provide the tools for using LANDSAT data for landcover 
classification. As identifed in the rest of this report, there have certainly been many painful 
lessons for project personnel about the limitations still present. But we have found ways to 
overcome these problems, and have created a landcover map of the Big Island that is the envy of 
the rest of the GIS community The final maps have generated great interest among the GIS 
community both in the maps themselves (everyone wants a copy), in acquiring LANDSAT and 
SPOT satellite imagery (SPOT has a cooperative project in the State of Hawaii), and in learning 
more about remote sensing. 
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APPENDIX A - Landsat TM Scenes and NOAA Data 


1) Scene Image: TM6346_90 

2) Scene Data Sheet: TM6346_90 

3) Scene Image: TM 6246_91 

4) Scene Data Sheet: TM6246_91 

5) Scene Image: TM6247_91 

6) Scene Data Sheet: TM6247_91 

7) Scene Image: TM6347_89 

8) Scene Data Sheet:TM6347_89 

9) Scene Image: TM6545_92 

10) Scene Data Sheet: TM6545_92 

11) NOAA Station Map 

12) NOAA Temperature Chart - June 1989 

13) NOAA Precipitation Chart - June 1989 

14) NOAA Temperature Chart - June 1990 

15) NOAA Precipitation Chart - June 1990 

16) NOAA Temperature Chart - October 1991 

17) NOAA Precipitation Chart - October 1991 
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Scene Image: TM6346 90 
Maui, NW Big Island 
June 1, 1990 
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EARTH OBSERVATION SATELLITE COMPANY 
*** Shipping Notice *** 



Ship Date: 08-feb-1994 

Customer Account #: 3000874 Customer Order #: 300087494010 

Invoice If : 94020050 Customer Ref ff : 0119312140063 

SHIP TO: 

V USGS/EROS DATA CENTER 

ATTEN: DATA MANAGEMENT/CEES 
MUNDT FEDERAL BUILDING 
SIOUX FALLS, SD, USA 

57198-000 

ATTN: LANDS AT DATA USER 

Phone # : 605-594-6511 



BILL TO: 

USDI/USGS 
EROS DATA CENTER 
CUSTOMER SERVICES 
SIOUX FALLS, SD, USA 

57198 

Product Product / . , Scene ID 

Order Number Code L# Path /Row /%Shift Acq-Date Mode Mux 


94010004-01 GLOBLT 4 j 063 j 046 j 00 90152 ( 01-jun-90 ) 1 0 

Maximum Full Scene Cloud Cover: 0 Bands Processed: 1234567 

Resample Method/Interval: CC /28.50 
Earth Ellipsoid: INTERNATIONAL 1909 

Map Projection: SPACE OBLIQUE MERCATOR " ^ 


Scene Data Sheet TM634O90 
Maui, NVV Big Island 
June 1, 1990 
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Scene Image: TM6246 91 
NE Big Island 
October 11, 1991 
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EARTH OBSERVATION SATELLITE COMPANY 
*** Shipping Notice *** 


Shi 

Customer Account #: 3000874 Customer Or 

Invoice #: 94020091 Customer Ref #: 119312150007 

SHIP TO: 

USGS/EROS DATA CENTER 
ATTEN: DATA MANAGEMENT/CEES 

MUNDT FEDERAL BUILDING 
SIOUX FALLS, SD, USA 

57198-000 

ATTN: LANDSAT DATA USER 

Phone #: 605-594-6511 


BILL TO: 

USDI/USGS 
EROS DATA CENTER 
CUSTOMER SERVICES 
SIOUX FALLS, SD, USA 

57198 


Product 
Order Number 

Product 

Code 

L# 

jpathl 

Row 

. Scene ID 

]%Shift Acq-Date Mode Mux 

94032003-01 

GLOBLT 

s| 

062 

046 

1 

00 91284 ( ll-oct-91 ) 1 0 

Maximum Full Scene 

Cloud ( 

1 

love r : 

/ 
3 C 

i Bands Processed: 1234567 


Resample Method/Interval: CC /28.50 



p Date: 15-feb-1994 
der #: 300087494032 


Earth Ellipsoid: INTERNATIONAL 1909 

Map Projection: SPACE OBLIQUE MERCATOR 


( 1 ) 


Scene Data Sheet TM6240 9I 
NE Big Island 
October II, 1991 


Appendix A r 




Scene Image: TM6247 91 
SE Big Island 
October II, 1991 
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EARTH OBSERVATION SATELLITE COMPANY 
*** Shipping Notice *** 


Ship Date: 16-feb-1994 

Customer Account #: 3000874 Customer Order #: 300087494032 

Invoice #: 94020104 Customer Ref #: 0119312150008 

SHIP TO: 

USGS/EROS DATA CENTER 
ATTEN: DATA MANAGEMENT/CEES 
MUNDT FEDERAL BUILDING 
SIOUX FALLS, SD, USA 

57198-0001 

ATTN: LANDSAT DATA USER 

Phone # : 605-594-6511 


BILL TO: 

USDI/USGS 
EROS DATA CENTER 
CUSTOMER SERVICES 
SIOUX FALLS, SD, USA 

57198 


Product 
Order Number 

Product 

Code 

L# 

^ath 

1 Row 

j Scene ID 

|%Shift Acq-Date 

Mode 

Mux 

94032004-01 

GLOBLT 

1 

■ 

062 

1 

047 

r _ ■ ~ 

00 91284 ( ll-oct-91 ) 

1 

0 


Maximum Full Scene Cloud Cover: 30 Bands Processed: 1234567 


Resample Method/Interval: CC /28.50 
Earth Ellipsoid: INTERNATIONAL 1909 

Map Projection: SPACE OBLIQUE MERCATOR W. 


( 1 ) 

Scene Data Sheet TM6247 91 
SE Big Island 
October II, 1991 
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EARTH OBSERVATION SATELLITE COMPANY 
*** Shipping Notice *** 


Customer Account #: 3000874 


Ship Date: 09-feb-1994 
Customer Order #: 300087494010 


Invoice #: 94020059 Customer Ref #: 0119312140061 


SHIP TO: 

USGS/EROS DATA CENTER 
ATTEN: DATA MANAGEMENT/CEE S 
MUNDT FEDERAL BUILDING 



CUSTOMER SERVICES 
SIOUX FALLS, SD, USA 

57198 


Product Product 

Order Number Code 


L# 


Scene ID 

Path /Row ftShift Acq-Date 


Mode Mux 


94010002-01 GLOBLT 


063 047 


Maximum Full Scene Cloud Cover: 


00 89165 ( 14- jun-89 ) 1 0 

Bands Processed: 1234567 


Resample Method/Interval: CC /28.50 
Earth Ellipsoid: INTERNATIONAL 1909 

Map Projection: SPACE OBLIQUE MERCATOR ^ 


( 1 ) 


Scene Data Sheet: TM6347_89 
S Big Island 
June 14, 1989 
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Scene Image: TM6545_92 
E Kauai, W Oahu 
April 1, 1992 
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EARTH OBSERVATION SATELLITE COMPANY 
*** Shipping Notice *** 

Ship Date: 16-feb-1994 

Customer Account #: 3000874 Customer Order # : 300087494032 

Invoice #: 94020100 Customer Ref #: 0119312150012 
SHIP TO: 

USGS/EROS DATA CENTER 
ATTEN: DATA MANAGEMENT/CEES 
MUNDT FEDERAL BUILDING 
SIOUX FALLS, SD, USA 

57198-0001 

ATTN: LANDSAT DATA USER 

Phone # : 605-594-6511 


BILL TO: 

USDI/USGS 
EROS DATA CENTER 
CUSTOMER SERVICES 
SIOUX FALLS, SD, USA 

57198 


Product 
Order Number 

Product 

Code 

L#J 

r 

Path 

Rowj 

Scene ID 

%Shift Acq-Date 

Mode 

Mux 

94032007-01 

GLOBLT 

4 1 

^065 

045 

00 92092 ( 01-apr-92 ) 

1 

0 


Maximum Full Scene Cloud Cover: 30 Bands Processed: 1234567 


Resample Method/Interval: CC /28.50 
Earth Ellipsoid: INTERNATIONAL 1909 

Map Projection: SPACE OBLIQUE MERCATOR ^ 




( 1 ) 


Scene Data Sheet: TM6545_92 
E Kauai, W Oahu 
April 1, 1992 
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APPENDIX B - Ancillary Field Data 


1) Sample Field Data Sheet 

2) GPS Sites Data Sheets 

3) Sample Cane Field Map and Cane Planting Data 

4) Hamakua Test Site Area 

5) Gorse in Pastures / Ohia-Koa Forest; Bands 4-5-3 

6) Gorse-ohia Classes over Raw Data 

7) Fountain Grass and Grazing Areas 

8) Kilauea Caldera, Volcanoes National Park 

9) Island of Hawaii 7.5 Minute Quad Index 

10) Soils Series Data from NRCS 



Jul 16 21:40 1995 ospwsOl : /usr/tmp/paaa26891 Page 1 


FIELD DATA VERIFICATION FORM - ISLAND OF HAWAII 

LANDSAT TM VEGETATION MAP 
STATE OF HAWAII -NASA HEADQUARTERS 

LANDSAT TM SCENES: 

HAWNW 63/46 JUNE1, 1990 HAWNE 62/46 OCTOBER 11,1991 

SOUTH POINT 63/47 JUNE14, 1989 KALAPANA 62/46 OCTOBER 11,1991 

SITE NUMBER ON MAP: DATE: 

ROAD MILE MARKER/TURNOFF: 

QUAD NUMBER: QUAD NAME: 

GPS NUMBER: GPS COORDINATES: N W 

ELEVATION: SLOPE/DIRECTION: 

LANDCOVER CLASS: 

ECOZONE DATA: 

COASTAL - DRY/MESIC/WET LOWLAND - DRY/MESIC/WET 

LOWLAND TO MONTANE - DRY/MESIC/WET MONTANE - DRY/MESIC/WET 
SUBALPINE - DRY/MESIC/WET ALPINE - DRY/MESIC/WET 

DOMINANT OR CO-DOMINANT GROUNDCOVER: 


MID-STORY GROUNDCOVER: 

UNDERSTORY GROUNDCOVER: 

GROUND COVER ON SOIL OR LAVA: 

TOPOGRAPHY : 

SOIL: TYPE/DEPTH LAVA FLOW: DATE/TYPE/NEW/WEATHERED/BROKEN/CRACKS/ROCKS 


NOTES /SKETCH/ PHOTOGRAPH ROLL -FRAMES : 


Sample Field Data Sheet 
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45 

A' a Dense Lichen 


4212400A 

E 

6 
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E 

17 
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4080321B 

I 

5 
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w 

38 
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u 
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C 

55 
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23 
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L 

54 
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412152 3A 

M 

31 
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Range 
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aa 

30 
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aa 
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P 

52 
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P 

7 
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P 
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N 
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N 
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H 
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Ohia-Dense 


4080220B 
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56 
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13 
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4080320A 

F 

15 

Ohia-Koa 


4080320C 
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43 
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4121323B 
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H 
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53 

Ohia-Olapa-Tree fern- Lava 
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Sorted by Land Cover 
Page 1 of 2 
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GPS Sites Data Sheet 
Sorted by Land Cover 
Page 2 of 2 
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1 Ohia-Lava 

2 Ohia-Dense 

3 Ohia-Lava 

4 Pahoehoe Lava 

5 Barren 
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7 Grass-Shrub Mamane-Naio 

8 Mamane-Mix Grass-Shrub 

9 Mamane-Mix Grass Shrub 
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11 Grassland-Range 

12 Koa-Ohia 
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Kiawe 


4080422A 

35 

Kiawe 


4080422B 

36 

Kiawe 


4080423A 

37 

Beach 


4080500A 

51 

Subalp Pukiawe-Ohelo-Aali 


4080S00A 

38 

Buf f elgrass-Kiawe 


4080500B 

39 

Grass land- Foun tain Grass 


4080501A 

40 

Grassland-Fountain Grass 

l 

4080501B 

41 

Grass land- Fountain Grass 


4080S02A 

42 

Ohia-Uluhe Fern 


4121323A 

43 

Ohia-Lava 


4121323B 

44 

Ohia-Uluhe-Tree Fern 


4121323C 

46 

Pine Forest-Grass-Fem 

Plantation 

4121401A 

47 

Ohia-A'a Lava 


4 12 14 2 0A 

48 

Shrub- Grass land 


4 121 42 1A 

49 

Grass land-Gorse- Lava 


4121421B 

50 

C i nder - Barr en 


4121423A 

52 

Grass -MaunaLani 


4121520A 

53 

Grass -Cactus - Shrub 

Rangeland 

4121521A 

54 

Eucalyptus - Plantation 

Forest Plantati 

4121S23A 

45 

A' a Dense Lichen 


421240 0 A 

54 

Lasuanna -Ohia- r em- Lava 


5031823A 

55 

Coastal Shrub-Grass-Lava 


5031900A 


GPS Sites Data Sheet 
Sorted by GPS ID# 
Page 1 of 2 
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GPS Sites Data Sheet 
Sorted by GPS ID# 
Page 2 of 2 
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Non-GPS Sites Data Sheet 
Sorted by Land Cover 
Page 1 of 1 
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Kilauea Caldera 
Volcanoes N.P., Big Island 
Oliia and M Faya 
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Volcano House, Kilauea Volocano: Bands 4, 2 and 1 ( r, g, b); 1 m/px 
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United States 
Department of 
Agriculture 


Natural 

Resources 

Conservation 

Service 


P. O. Box 636 
Kealakekua, HI 
96750-0636 
Phone (808) 322-2482 
Fax (808) 322-3735 


General Soils Map — Island of Hawaii 
Soil Series / Map Units 
Draft 6 - 5/95 
file: Gen-map-06 

Robert Gavenda, Soil Survey Project Leader 
Note: 

1) The proposed classifications presented here are best estimates and are based on limited 
(often extrapolated) lab data. The classification presented here assumes that Hydrudands will 
key out before Melanudands, and Folist great groups will be based on moisture regimes. 

2) All new series and changes in classification in established series are in italics. 


Map Units 

1 - Lava flows, a' a and pahoehoe undifferentiated 
Published series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Rock land 

Proposed series/ miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Cinder land 
Typic Ustifolists 
Lithic Ustifolists 
Typic Udifolists 
Lithic Udifolists 

2 - Lava flows - Cinder land 
Published series/ miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Cinder land 

Proposed series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Cinder land 

Nenenui: pumiceous, isomesic Lithic Ustorthents 
?????: euic, isomesic Typic Ustifolists 

?????: euic, isomesic Lithic Ustifolists Soils Series Data from NRCS 

3 - Lava flows - Haplustolls - Torriorthents (aridic, hyperthermic) Appendix B #10 



Palapalai: medial, isothermic Humic Haplustands 

11 - Fulvudands (udic, isothermic) 

Published series/miscellaneous land types: 

Niulii: thixotropic, isothermic Hydric Dystrandepts 
Proposed series/ miscellaneous land types: 

Niulii: hydrous , isothermic Hydric Fulvudands 

12 - Andaquepts - Endoaquands (udic, isothermic) 

Published series/miscellaneous land types: 

Kehena: thixotropic, acid, isothermic Aerie Andaquepts 
Amalu: fine, mixed, acid, isomesic Histic Placaquepts 
Tropaquepts 
Rough broken land 

Proposed series/miscellaneous land types: 

Kehena: hydrous , acid, isothermic Aerie Andaquepts 
?????: hydrous, isothermic Histic Endoaquands 

13 - Hydrudands - Placudands - Placaquepts (perudic, isomesic) 
Published series/miscellaneous land types: 

Manahaa: thixotropic, isomesic Hydric Dystrandepts 
Kahua: thixotropic, isomesic Typic Placandepts 
Maile: thixotropic, isomesic Hydric Dystrandepts 
Amalu: fine, mixed, acid, isomesic Histic Placaquepts 
Proposed series/miscellaneous land types: 

Kahua: hydrous, isomesic Hydric Placudands 
Maile: hydrous, isomesic Typic Hydrudands 
Amalu: fine, mixed, acid, isomesic Histic Placaquepts 

14 - Haplustolls (aridic, thermic) 

Published series/miscellaneous land types: 

Puu Pa: medial-skeletal, isothermic Ustollic Eutrandepts 
Kamakoa: ashy over loamy isothermic Mollic Ustifluvents 
Proposed series/ miscellaneous land types: ajmc. 

Puuhinai: medial-skeletal, thermic Torrvrtia Haplustolls 
?????: medial, thermic 'Ferrrrrtic Haplustolls 
Kamakoa: ashy over loamy, thermic Mollic Ustifluvents 
Popoo: medial-skeletal, thermic Typic Haplustolls 
Lahuipuaa: sandy-skeletal, thermic Typic Torriorthents 

15 - Haplustands (ustic, isothermic) 

Published series/ miscellaneous land types: 

Puu Pa: medial-skeletal, isothermic Ustollic Eutrandepts 
Waikaloa: ashy, isothermic Ustollic Eutrandepts 
Waimea: medial, isothermic Typic Eutrandepts 
Kikoni: medial, isothermic Typic Eutrandepts 
Kilohana: ashy, isomesic Typic Ustipsamments 
Proposed series/ miscellaneous land types: 

Waikaloa: medial, isothermic Typic Haplustands 
Waawaa: medial, isothermic Typic Haplustands 
Waimea: medial, isothermic Typic Haplustands 
Kikoni: medial, isothermic Typic Haplustands 

16 - Ustifolists - Haplustands - Ustorthents (ustic, isothermic) 
Published series/miscellaneous land types: 

Lava flows, a' a 
Lava flows, pahoehoe 
Very stony land 



Hanipoe: medial, isomesic Thaptic Haplustands 
Laumaia: medial, isomesic Humic Haplustands 
Cinder land 
Lava flows, a' a 

21 - Fulvudands (udic, isomesic) 

Established series/miscellaneous land types: 

Umikoa: medial, isomesic Entic Dystrandepts 
Hanipoe: medial, isomesic Typic Dystrandepts 
Puu Oo: thixotropic, isomesic Hydric Dystrandepts 
Puukala: thixotropic-skeletal, isomesic Hydric Lithic Dystrandepts 
Proposed series/ miscellaneous land types: 

Umikoa: medial, isomesic Hydric Fulvudands 

Hanipoe: medial, isomesic Hydric Fulvudands 

Puu Oo: hydrous , isomesic Hydric Fulvudands 

?????: hydrous-skeletal, isomesic Hydric Lithic Fulvudands 

22 - Hydrudands (perudic, isomesic) 

Established series/miscellaneous land types: 

Piihonua: thixotropic, isomesic Typic Hydrandepts 
Akaka: thixotropic, isomesic Typic Hydrandepts 
Proposed series/ miscellaneous land types: 

Piihonua: hydrous, isomesic Typic Hydrudands 
?????: dysic, isomesic Typic Troposaprists 
Piha: ????? 

23 - Hydrudands (perudic, isothermic) 

Established series/ miscellaneous land types: 

Kaiwiki: thixotropic, isothermic Typic Hydrandepts 
Akaka: thixotropic, isomesic Typic Hydrandepts 
Keei: euic, isothermic Lithic Tropofolists 
Proposed series/miscellaneous land types: 

Kaiwiki: hydrous, isothermic Acrudoxic Hydrudands 
Akaka: hydrous, isothermic Acrudoxic Hydrudands 
?????: dysic, isothermic Typic Troposaprists 

24 - Hydrudands (perudic, isohyperthermic) 

Established series/miscellaneous land types: 

Hilo: thixotropic, isohyperthermic Typic Hydrandepts 
Proposed series/miscellaneous land types: 

Hilo: hydrous, isohyperthermic Acrudoxic Hydrudands 

25 - Lava flows - Ustifolists - Ustorthents (ustic, isothermic) 
Established series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Cinder land 

Kekake: dysic, isomesic Lithic Tropofolists 
Mawae: euic, isomesic Typic Tropofolists 
Kilauea: ashy-skeletal, isothermic Typic Ustorthents 
Kona: euic, isothermic Lithic Tropofolists 
Proposed series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Cinder land 

Kilauea: ashy-skeletal, isothermic Typic Ustorthents 
?????: euic, isothermic Lithic Ustifolists 

26 - Lava flows - Udifolists (udic, isomesic) 

Established series/miscellaneous land types: 



Olaa: thixotropic over fragmental, isohyperthermic Typic Hydrandepts 
Panaewa: thixotropic-skeletal, isohyperthermic Lithic Hydrandepts 
Proposed series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 

Papai: euic, isohyperthermic Typic Udifolists 

Keaukaha: euic, isohyperthermic Lithic Udifolists 

Opihikao: euic, isohyperthermic Lithic Udifolists 

Malama: euic, isohyperthermic Typic Udifolists 

Olaa: hydrous-skeletal , isohyperthermic Typic Hydrudands 

Panaewa: hydrous-skeletal, isohyperthermic Lithic Hydrudands 

30 - Perudifolists - Lava flows (perudic, isothermic) 

Established series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Keei: euic, isothermic Lithic Tropofolists 
Kiloa: euic, isothermic Typic Tropofolists 
Proposed series/miscellaneous land types: 

Lava flows, a'a 

Lava flows, pahoehoe 

?????: dysic, isothermic Lithic Perudifolists 

?????: dysic, isothermic Typic Perudifolists 

31 - Perudifolists - Hydrudands - Lava flows (perudic, isomesic) 
Established series/ miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Keei: euic, isothermic Lithic Tropofolists 
Kiloa: euic, isothermic Typic Tropofolists 
Lalaau: euic, isomesic Typic Tropofolists 
Kahaluu: euic, isomesic Lithic Tropofolists 
Puukala: thixotropic-skeletal, isomesic Hydric Lithic Dystrandepts 
Piihonua: thixotropic, isomesic Typic Hydrandepts 
Proposed series/ miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 

Kahaluu: dysic, isomesic Lithic Perudifolists 
Lalaau: dysic, isomesic Typic Perudifolists 
?????: hydrous-skeletal, isomesic Lithic Hydrudands 
Piihonua: hydrous, isomesic Typic Hydrudands 

32 - Ustorthents - Haplustands - Lava flows (ustic, isothermic) 

Established series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 

Kilauea: ashy-skeletal, isothermic Typic Ustorthents 
Heake: medial, isothermic Lithic Dystrandepts 
Proposed series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 

Kilauea: ashy-skeletal, isothermic Typic Ustorthents 
Heake: medial, isothermic Lithic Haplustands 

33 - Ustifolists - Haplustands - Lava flows (ustic, isomesic) 

Established series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 

Hanipoe: medial, isomesic Typic Dystrandepts 
KekaJke: euic, isomesic Lithic Tropofolists 



Naalehu: medial, isohyperthermic Humic Haplustands 

38 - Haplustands (ustic, isothermic) 

Established series/ miscellaneous land types: 

Kamaoa: medial, isothermic Typic Eutrandepts 
Proposed series/miscellaneous land types: 

Kamaoa: medial, isothermic Typic Haplustands 
Keaa: medial-skeletal, isothermic Lithic Haplustands 

39 - Haplustands (ustic, isohyperthermic) 

Established series/miscellaneous land types: 

Pakini: medial, isohyperthermic Entic Eutrandepts 
Kaalualu: medial-skeletal, isohyperthermic Ustollic Eutrandepts 
Naalehu: medial, isohyperthermic Typic Eutrandepts 
Proposed series/miscellaneous land types: 

Pakini: medial, isohyperthermic Typic Haplustands 
Kaalualu: medial-skeletal, isohyperthermic Humic Haplustands 
Naalehu: medial, isohyperthermic Typic Haplustands 

40 - Lava flows - Ustifolists (ustic, isohyperthermic) 

Established series/ miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 

Punaluu: euic, isohyperthermic Lithic Tropofolists 
Kaimu: euic, isohyperthermic Typic Tropofolists 
Kainaliu: medial-skeletal, isohyperthermic Typic Eutrandepts 
Proposed series/miscellaneous land types: 

Lava flows, a'£ . 

Lava flows, pahoehoe 

Punaluu: euic, isohyperthermic Lithic Ustifolists 

Kaimu: euic, isohyperthermic Typic Ustifolists 

Kainaliu: medial-skeletal, isohyperthermic Typic Haplustands 

41 - Ustifolists - Haplustands - Lava flows (ustic, isohyperthermic) 
Established series/miscellaneous land types: 

Lava flows, a’a 
Lava flows, pahoehoe 

Punaluu: euic, isohyperthermic Lithic Tropofolists 
Kaimu: euic, isohyperthermic Typic Tropofolists 
Kainaliu: medial-skeletal, isohyperthermic Typic Eutrandepts 
Waiaha: medial-skeletal, isohyperthermic Lithic Eutrandepts 
Proposed series/ miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 

Punaluu: euic, isohyperthermic Lithic Ustifolists 
Kaimu: euic, isohyperthermic Typic Ustifolists 
Kainaliu: medial-skeletal, isohyperthermic Typic Haplustands 
Waiaha: medial-skeletal, isohyperthermic Lithic Haplustands 

42 - Ustifolists - Lava flows (ustic, isothermic) 

Established series/miscellaneous land types: 

Lava flows, a'a 
Lava flows, pahoehoe 
Puna: euic, isothermic Typic Tropofolists 
Kona: euic, isothermic Lithic Tropofolists 
Proposed series/miscellaneous land types: 

Lava flows, a'a 

Lava flows, pahoehoe 

Kapua: euic, isothermic Typic Ustifolists 



Hokukano: hydrous, isomesic Lithic Haplustands 
Kealoha: medial-skeletal, isomesic Typic Haplustands 
Hapuu : loamy-skeletal, isomesic Ustic Humitropepts 
Nenenui: pumiceous, isomesic Lithic Ustorthents 
Ohianui: pumiceous, isomesic Lithic Ustorthents 



APPENDIX C - Classification System 


1) Raw Data with Jacobiveg Level II Overlay 

2) Sample Class Name List and Map 

3) Hilo Urban Areas with Roads 

4) Signature File and Legend 

5) Island of Hawaii - Landsat TM Subtropical Hierarchical Classification System 



Raw Data with Jacobi Veg Overlay 
Saddle Road, Big Island 
Bands 7-5-3, Jacobi Level II 
GPS sites in yellow text 

Appendix C #1 








































Legend 


Signature File 

L tend 


Unclassified 

Unclassified 


Dense Wet Ohia-Koa vv Hapu'u 
Tree fern. Shrub, Grass; HVNP 

Lowland to Montaine 
Dense Ohia Koa Wet 
Forest 


Recent Lava, Mauna Loa Summit 

Bare Lava 


Lowland Wet Rain Forest (Mauulio 
Stream Gulch) 

Lowland Dense Mixed 
Wet Forest 


Montane Ohia-Wcather Lava 1855- u 
Mx Gras, Mat Fern GPS408022 1 B 

Montaine Open Ohia 
Mesic Forest 


Bare Soil - Hamakua Sugar Co Field 
20470 W Kaumoali Stream 

Bare Soil 


Older Lava, Mauna Loa 

Bare Lava 


Recent Lava; 1843 1935 flows- W 
Pu'u Nene; GPS4080223A; A’A Lava 

Bare Lava 


Recent Pahoehoe Lava 
GPS4080222A 

Bare Lava 


Shadow 

Shadow 


Water (hawnw_sc!) 

Water 


Mature Sugar Cane- Hamakua Sugar 
Co Field 54060/51310 E Paauilo 

Sugar Cane 


Young Sugar Cane-Soil- Hamakua 
Sugar Co Field 43200; N Paauilo 

Sugar Cane 


Bare Soil, Hamakua Suear Field 
22300/203 10; E Honokaa on H-19 

Bare Soil 


Dry Koa Community-w Exotic 
Grassland; Pasture GPS4080400F 

Montaine Open Mesic 
Koa Forest (Grazed) 


Bare Lt. Soil- Hamakua Fields 
24390/24630/3 1370; SE of Honokaa 

Bare Soil 


Coastal Kiawe Forest GPS4080423A 

Coastal Open Dry Kiawe 
Forest 


Montane Dense Mesic Ohia-Uluhe 
Fcm, GPS4121323A 

Montaine Dense Mesic 
Ohia Forest 


Bare Soil Hamakua Field 88600; SE 
Lapahoehoe W Maauilo Stream 

Bare Soil 
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Montane Dry Koa Community- w X 
Grass-Mx Shrub; Rang 
GPS4080400E 


Barren-Cinder Mauna Kea 


Cinder- E. Mauna Kea Summit by 
Pu'u Wekiu 


Wet Koa-Ohia Forest vv/ other Native 
Trees, Shrub, Trcefem 


Barren-Cinder 1, Mauna Kea; 
GPS4121423A 


Beach GPS4080500A 


Subaipine Dry Na Shrub-Lava; 
Pukiawe, Oheio, A'ahi 
GPS412I500A 


Lovvlan Rain Forest-w Ohia-Koa, 
Hapu'u Treefm, Gras GPS4080320A 


Low-Montane Rain Forest- vv Koa, 
Hapu'u <80% Canopy 
GPS4080320B 


Montane Wet Ohia-Koa >65% 
Canopy GPS4080320C 


Montane Dry X Tropic Ash 
Plantation 80% Canopy 
GPS408032 1 A; Kal 


Montane Dry X Tropicl Ash 
Plantati-B Poka Liana 
GPS4080321B;Kal 


Montane Dry Koa-w B. Poka Liana 
on trunks-canopv GPS4080322A 


Montaine Open Dry Koa 
Forest (Grazed) 


Bare Cinder 


Bare Cinder 


Montame Open Ohia Wet 
Forest 


Bare Cinder 


Bare Sand Beach 


Subaipine Open Pukiawe 
Oheio Aalii Shrubtand 


Lowland Dense Ohia Koa 
Wet Forest 


Lowland to Montaine 
Dense Ohia Koa Wet 
Forest 


Montaine Open Ohia Wet 
Forest 


Forest Plantation 
fFraxinus) 


Montane Dry Koa-with Banana Poka 
Liana GPS4080322B 


Coastline 


Montane Mesic Open Koa-w B. P. 
Liana, Fun-Grass GPS40S0322C 


Forest Plantation 
(Fraxinusl with Banana 


Poka 


Montame Open Koa 
Mesic Forest with Banana 
Poka 


Montaine Open Koa 
Mesic Forest with Banana 
Poka 


Coastal Surf 


Montaine Open Koa 
Mesic Forest with Banana 
Poka 
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Montane Dry Koa-w Banana Poka 
Liana-Grass GPS4080400A 

Montame Open Koa 
Mesic Forest with Banana 
Poka 


Wet Ohia Forest- w Treefem-Native 
Shrub undersl E-GPS4080220B 

Montaine Open Ohia Wet 
Forest 


X Gras-Mx Shrub w Fountain 
Grass, Cacti, Mullein; SE 
GPS4121523A 

Lowland Dry Open 
MixedShrub lands with 
Fountain Grass 


Gorse-Mixed Grasses, Soil NW 
GPS4080402C 

Montaine Wet Grasslands 
with Gorse 


Montane Wet Exotic 
Grass-Rangeland GPS4080402A 

Montaine Wet Grasslands 
(grazed) 



Sugar Cane, Hamakua Sugar Co. 
Field 20330; SE of flonokaa 

Montane Wet X Gorse-Nat 
Mamane-X Grass, Cinder, h ne 
GPS4080402C 


Abandoned Sugar Cane, Hamakua 
Field 11030, EofWCpio 

Montane Dry X Grass-Mullein- w 
Koa-Mamane; Rangeian 
GPS4080403A 

Subalp Dry Mamane Woodland-Mix 
Grass, Shrub, Cindr GPS4080301A 


Sugar Cane 


Montaine Wet Grasslands 
with Gorse 


Sugar Cane 


Montaine Dry Open Koa 
Mamane Forest (grazed) 


Montaine Dry Open 
Koa Mamane Forest 



Coastal-Lowland Exotic Lowland Dry Open 

Grassland-Lava with Fountain Grass Grasslands with Fountain 

Grass 


Cloud Shadow (forest Wof 
Oncma-River) 

Montane Dry Mamane Woodlan-w 
Mx Grass-Shrub 50% Can 
GPS4080300B 


Montane Dry Mamane Woodlan-w 
M\ Grass-Shrub 50% Can 
GPS4080300C 

Lowlan-Montane Mamane-Naio 
Woodland-vv Mix Shrub-Grass 


Shadow 


Montaine Dry Open 
Mamane Forest 


Montaine Dry Open Forest 


Montaine Dry 

OpenA famc)ie/?\ aio Forest 



Subal D Gras-S^rub-Cinr-w Subalpine Dry’ Open 

Desch,K. Hade, Pul:, Aali Mixed Shrub land 

GPS4080300A 
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Low land Dry Grassland w Fountain 
Grass GPS4080501 A 

Lowland Dry Open 
Grasslands with Fountain 
Grass 


Barren-Cinder 2, Mauna Kea, 
GPS4121423A 

Cinder 


Barren Cinder 3, Mauna Kea; 
GPS4121423A 

Cinder 


Dense Young Ohia-Uluhe fern- w 
dead tree trunks, GPS4 12 1 323 A 

Montaine Wet Dense Ohia 
Forest with Uluhe 


Montane Nativ Shrubs-Lava; 
Pukiawe, Ohelo, A'ali'i 
GPS412 1500A 

Montaine Mesic Open 
Mixed Native Shrubland 


Dr}' Gras-Shrub w Fountain 
Gras, Cacti -r/or Mullein E 
GPS4I21523A 

Lowland Dry Open 
Grasslands with Fountain 
Grass 


Coastal Kiawe Forest >70% Canopv 
GPS4080422B 

Coastal Dry Dense Kiawe 
Forest 


Dry Grass-Rock, Shrub, P. Pear 
Cactus, Mullein SE GPS4 1 2 1 523 A 

Lowland Dry Open 
Grasslands 


Coastal Mix Grass- A' A Rock w open 
dead Kiawe Forest GP S4080500B 

Coastal Dry Open Kiawe 
Forest 


Lowland Dry Exotic Eucalyptus 
Forest Plantation GPS4 1 2 1*523 A 

Forest Plantation 
fEucalvptus') 


X Grass-Shrub, P.P. Cactus, Mullein, 
Range NE Saddle Rd./H- 1 90 

MontaineDry Open 
Shrub lands 


Lava Flow 1 899- w Montane Nat 
Shrubs Pukiawe-Ohclo 
GPS4121421 A 

Montaine Mesic Open 
Shrub lands 


Montane Grass lan- w Molasses-B 
Beard-Broomsedge-Pili Steam Vent 

Montaine Mesic 
Grasslands 


Class 34-CIoud-Urban from 
Unsupervised 40 Classification 

Cloud 


Montane Dense Wet 
Ohia-Shrub-Ama'u Fem-L Flow 1 855 
GPS4080220A 

Montaine Wet Dense 
O/wa/Shrub Forest on 
Lava 


Class 36 -Cloud from Unsup 40 
Classification 

Cloud 


Class 37-Cloud from Unsup 40 
Classification 

Cloud 
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Montane Wet Ohia- Weathered Lava 
Flow of 1881; GPS4 1 2 1 323B 


Montane Dense Wet Ohia-Hapu'u 
Tree fem-Uluhe fern; GPS4 1 2 1 323C 


Montane AW Lava Flow 1 984-Dense 
Lichen GPS4121400A 


Montane Mcsic Open Ohia-Ululie 
Fern (6 in.-2 f ); GPS4I21420A 


Montane Slash or Loblolly 
Pine-Forest Plantation; 
GPS4121401 A 


Coastal Grassland-Molasses-B. Beard 
w Pluchea-Guava-Rose Apple 


Cloud shadow-Kalapana 


Montaine Wet Dense Ohia 
Forest on Lava 


Montaine Wet Dense 
Ohia^Hapuu Forest 


Montaine Wet Dense 
Lichen on Lava 


Montaine Mesic Open 
Ohia/Uhihe Forest 


Forest Plantation (Pine) 


Coastal Mesic Mixed 
Slxrub/ Grasslands 


Cloud Shadow 



Cloud Shadow S of Keanakolu Rd. Cloud Shadow 


Montane X Grassland-Mullein- w 
Koa SW GPS4080403A MK Rd 


Macadamia Orchard, Mauna Loa 
Gardens 


lava 5 kileauea crater-Kalapana 
(UCL40) 


Roseapple ( Vi ve)-Guava- Strawberry 
Guava w Molasses-Mix Grass 


Cinder- Ash; S W Kilauea Crater Rim 
Drive 


Montane Dry Koa-w Banana Poka 
Liana, Grass GPS4080400A; hawmv 


Coastal Wet Grassland- Shrub - 
Molasses w Pluchea; burn 83/6-91 


Water (hawneuc! merge 9-16) 


Montane Mesic t . gi Pine Forest 
Plantation GPS 4080401 A; hawmv 


Exotic Grass-Rangeland 
GPS4080403A; hawnw 


Coastal Kiawe Forest-Soil 
GP S40804 2 3 A;hawmv 


Montaine Open Mesic 
Koa Forest (grazed) 


Macadamia Orchard 


Bare Lava 


Coastal Mesic Mixed 
Shrub/Grasslands 


Cinder 


Montaine Open Koa Dry 
Forest with Banana Poka 


Coastal Wet Mixed 
Shrub/Grasslands 


Open Water 


Forest Plantation (Sugi) 


Montaine Mesic 
Grasslands (grazed) 


Coastal Dr}' Open Kiawe 
Forest 
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Montane Dry Koa Community- w Mx 
Gras-Shrub; Rang GPS4080400E;n\v 


Montane Dry Koa Community- w Mx 
Gras-Shrub Range GPS4080400F n\v 


Cloud from sup40; Spont 


Deep Water; Spomt 


Cloud Shadow; Spoint 


Cloud shadow over cane; Spoint 


Lowland Wet Koa-Ohia w Hapu’u 
Trefem-w Shrub-Gras-CHdemia-SG 


Montaine Dry Open Koa 
Forest (grazed) 


Montaine Dry Open Koa 
Forest (grazed) 



Shadow 



Hilo Bay 


Kealakekua water 


Dry Koa-Banana Poka Liana with 
Grass GPS4080400A 


Lowland Grassland-Mix Shrub; 
Fountain Grass, Koa Haole-Pluchea 



Melastomes 


Coastal Grassla-Mix Shrub-Lava; 
BufTel Grass-Pluchea-Koa Haole 


Montane Mesic vog Ohia-Uluhe-w 
BBeard Mol as Brooms Gras 


Halema'uma’u-Kalapana Volcanic 
Gases- Solfataras sulfer fumes 


Monta D Ohia Wodlnd-w Ohia, 
Aali-Ulei-Puk-Ohelo-Cindr; Klookout 


Montane Dry Shruhland-w A’ali'i, 
Pukiawe, Ohelo, Cinder; USGS M 


Shadow 


Lowland Wet Open 
Koa' Ohia Forest with 
Ilapuu and Clidemia 


Open Water 


Open Water 


Montaine Open Koa 
Mesic Forest with Banana 
Poka 


Lowland Dry Open 
Shrub/Grasslands with 
Fountain Grass 


Lowland Mesic Dense 
Shrublands with 
Melastomes 


Coastal Dry Open 
Shrub/Grasslands on Lava 


(either Montaine Mesic 
Open Ohia 
Forest/S hrub land OR 
Vog) 




Montaine Dense Ohia 
Forest/Shrubland on Lava 


Montaine Open shrub land 
on Lava 



Class 1-Kilauea Iki, halema'uma'u Bare Lava 


U ater-Rescrvoir S W Hilo 


Halema’uma’u-Kilnuea Iki 


Open Water 


Bare Lava 
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LLJaaaju) 

F dsL QfvSL 


/ n j . . _ 


Ks/IS Cj's4- 


|0l 0{iAi{ \o 3^0 

to SX& 



1 'T-L 

Pahochoe Lav a Flow - 1 974 

Bare Lava 


Lovvlan Dry Koa Community w 
Molasses-B. Beard Grass; Ililina Pali 

Lowland Dry Open Koa 
f* orest/Grassl and 


Kilauea Crater 

Bare Lava 


Lowland Dry Grassland- with 
Molasses Grass, Thatch, Pili 

Lowland Dry Grassland 


Mont dead Ohia 

Woodland- Shmb-Grass-Fem-Bcny; 
E Kipuka Puaulu 

Montaine Mesic Open 
Ohia Forest/Shrubland 


Class 2-1 Ialema’uma'u Kilauea, Puu 
O'o, Hilina Pali 
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Montaine Wet Dense 
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Firetree (Mynca faya)- Volcano House 

Montaine Wet Dense 
Mixed Forest with Firetree 


Macadamia Orchard_SPoint 
Scenefsmall for hnw) 

Macadamia Orchard 


Silk Oak-Grass_SPoint 

Montaine Mesic Dense 
Mixed Forest with Silk 
Oak 


Low-Mo Dry Ohia-Mfava Forcst-w 
Aali, Amau, fern, Cin-Ash 2,500' 

Lowland Dry Open Ohia 
Forest/Shrubland on Lava 


Lowl D vog Shrubl-w Aalii, few 
Ohia, Molas Grass; hairpin turn 

Lowland Dry Open 
Shrub/Grassland on Lava 

. • j 

Low D vog Ohia Forest- w Aaii, Mola 
Gras, Spars Ohia E hairtum 

Lowland Mesic Dense 
Ohia Forest/S hrub land 


Macadamia Orehard-C. Brewer, Kail 
Sugar Co. Pahala flield size) 

Macadamia Orchard 


Class 6-IIalema uma'u Kilauea, 
Mauna Lua Flow to coast 

Bare Lava 


lava 7-Kilauea Crater and 1 9?? flow 
j Kau Desert, Kona 19?? 

Bare Lava 
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.73 

!S. L0 T l8 ? d ^ ?’*" MiXCd Forest - Gr «sland * with Molasses-Razor Grass, 
Pluchea, lantana, Oh i a- Ha i a 

.82 

Lowland Wet Grassland - with dense Molasses Grass, Pluchea, Bamboo Orchid 
Lowland Uet Dense Koa-Ohia Forest - with Hapu'u Treefern, Mixed Shrub-Grass 

T iC °T e MiXed FOreSt • “ ith 0hia - Hata ' Roseapple, Wiliwili, 

Loilu' ' I!!T’ T 8 ' M ° nkey POd ' Slrawberry Guava - Guava . Java Plus; 

Ulei ' Pluchea ^ A, ^ r alian Palms, Christmasberry, Silk Oak, Hapu'u Treefern, 
me,, Pluchea, mixed sedges-grasses with Ohe-Uki, Maupaka, ferns Ama'u-Uluhe 

Lowland Uet Open Koa Forest - with Hapu'u Treefern, Banana Poka and other Lianas 


Montane Uet Ohia-Shrub - with Pukiawe, Ohelo, Orchid on weathered lava 1855 flow 

Montane Uet Open Ohia-ICoa Forest - with Hapu'u Treefern, mixed shrub grass 

Montane Uet Grassland - with Gorse, mixed grasses 

Montane Uet Grazed Grassland - with Kikuyu, mixed grasses 

Montane Uet Open Shrubland-Lava - with Ama'u, Ohia-Deschanpsia ecozone edge 

Montane Uet Lichen on 1984 Lava Flow - with Pioneer Ohia along road-ecozone edge 


.98 

Lowland Mesic Dense Shrubland - with Melastomes 

• 14 

Montane Mesic Open Grazed Koa Forest - with Kikuyu-mixed grasses, Rubus, trunks 

Montane Mesic Open Grazed Grassland - with Molasses Grass, Bush Beard Grass, 
PHi, Broomsedge, Bamboo Orchid on weathered Pahoehoe Lava 

Montane Mesic Open Ohia-Uluhe - with Bush Beard Grass, Bamboo Orchid 

• Ilf 

Montane Mesic Mixed Forest - with Fire Tree-Ohia, Molasses Grass, mixed grasses 
Montane Mesic Dense Mixed Forest - with Silk Oak, Ohia-Koa, Beard Grass, Ulei 
Sifcalpine Mesic Open Shrubland - with Pukiawe, Ohelo, Lepo-nene on lava 


.16 

Coastal Dry Kiawe Forest - with Kiawe, Buffelgrass. Pluchea on weathered lava 

Coastal Dry Open Grass-Shrubland - with Buffelgrass, Pluchea, Lantana, Koa hade 

Lowland Dry Open Shrub- Grass l and - with Koa haole, Pukiawe, Prickly Pear Cactus, 
Lantana, dense Fountain Grass 
.110 


Lowland Dry Dense Grassland - with Molasses, Pili, Bush Beard Grass, A'ali'i 

Lowland Dry Open Mixed Forest-Shrubland - with Molasses Grass 
.49 

Montane Dry Open Namane-Naio Uoodland and Forest - with scattered Ohia, Koa 

mixed grasses, mixed shrub, on weathered lava 

.43 

Montane Dry Open Koa-Mamane Forest - with A'ali'i, Pukiawe, Bush Beard Grass, 
Uki, Deschampsia 



.61 

Montane Dry Open Grazed Grass -Shrub Land - with Fountain Grass, mixed grasses. 

Mullein, Prickly Pear Cactus, Koa haole, l an tana, A'ali'i, Koa trunks 

.50 

Subalpine Dry Open Grass- Shrub Land - with Dechampsia, Pukiawe, Ohelo on lava 
.25 

Subalpine Dry Native Shrubland - with Pukiawe, Ohelo, A'ali'i on lava 
.44 

Subalpine Dry Open Mamane Forest - with Pukiawe, Ohelo, Deschampsia on lava 


.2 

Bare Lava 
.5 

Bare Soil 

.20 

Cinder Fields- Ash 
.9 

Shadow, Water 

.11 

Sugar Cane - data from Hamakua and Kau Sugar Companies 
.77 

Macadamia Orchard and Coffee 
.64 

Cloud, Snow, Volcanic Gases, Urban 



Island of Hawai'i-Landsat TM 

^ Subtropical Hierarchical Land Covor Classification 

Physical Ecosystems] | Biological Ecosyst^ Manmade Ecosysi^ 















APPENDIX D - Mosaic Process 


1) Landsat TM Scenes Before Mosaic 

2) Mosaic Model in Imagine 

3) Sample Area of Interest for Mosaic 

4) Raw Data Mosaic of the Island of Hawaii with Maui 

5) Classified Mosaic and Legend of the Island of Hawaii 
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APPENDIX E - LSVIEW Quick Tour 


1) Sample Screens from LSVIEW with Explanations 

1: LSVIEW Opening Full Screen Layout 

2: LSVIEW Main Menu 

3: LSVIEW Location Map Menu 

4: LSVIEW Detail of Location Map for Mauna Kea 

5: LSVIEW Full Screen Showing Mauna Kea Location 

6: LSVIEW Map Display Menu 

7: LSVIEW Map Display Menu with TM 7-5-3 Mosaic Active 

8: LSVIEW Full Screen with TM 7-5-3 Mosaic of Mauna Kea Area 

9) LSVIEW Main Menu to Show Site Photos 

10: LSVIEW Full Screen Showing Site Photos Near Mauna Kea 
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LSVIEW OPENING FULL SCREEN LAYOUT 



Custom View! User defines new window to zoom 

into on Detail Map. 
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LSVIEW MAIN MENU 
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LSVIEW LOCATION MAP MENU 
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LSVIEW QUICK TOUR - page 5 



Clicking ‘Map Display...’ in the Main Menu brings up 
this Map Display Menu. 
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LSVIEW MAP DISPLAY MENU 



TM 7-5-3 Mosaic has been made active so it will be 
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LSVIEW QUICK TOUR - page 7 







The last button to illustrate is ‘Show site photos!’ 
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LSVIEW MAIN MENU TO SHOW SITE PHOTOS 





LSVIEW Quick Tour 
Page 10 of 10 
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LSVIEW FULL SCREEN SHOWING SITE PHOTOS NEAR MAUNA ICEA 




APPENDIX F - Change Detection 


1) MSS 1977 Raw Data Mosaic 

2) Change Detection from the University of Hawaii 

1: Lava Flows at Hawaii Volcanoes National Park 

2: Gorse and Koa-Banana Poka 

3: Eucalyptus Plantation Keyhole 

4: Change in Sugar Cane Cultivation 

5: New Golf Courses 

6: Mauna Kea Forest Reserve Area 

3) List of MSS and TM Band Combinations Revealing Significant Change 

4) TM Raw Data Mosaic 



0 



MSS 1977 Raw Mosaic 
Big Island 
Acquired same day 
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Lava Flows 


Lava Flow of 
1970- 1974 





k,. , 

y-’ 



H * w Lava - 
Kalapana 
1983 to date 






The top image is the 1 977 MSS data. The bottom image is the 1991 TM data combined with the 
MSS data depicting the change in this land area. New lava flows are shown in green in die upper 
right portion of this image. Note: Landsat TM is currently the only satellite that provides shortwave 
bands 7 and 5 to detect active lava and provide maximum separation of lava, soil, and vegetation. 


Digital data indudes 4 scenes from Landsat TM (June 1989, 

June 1 990, and Oct. 1991) registered with 2 scenes from 
Landsal MSS (Oct. 1977). AH geographic databases of the 
State of Hawaii Geographic Information System are in UTM 
Zone 4, although (he Island of Hewai k i is in UTM Zone 6. This 
image is projected into UTM Zone 4 to comply with the State 
of Hewai 'i database design. Base is adapted from Hswai 'i Island 
USGS 1:24.000 topographic maps. Old Hawaiian Datum is based 
on Clarke 1866. MSS processing by Ed Petteys, Chris Hogan, 

Pat Costales, Mike Buck, and Craig Tasaka. TM processing by 
Chris Hogan, Ron Cannarefla, Leo Ikebara, and Roger Imoto. 


This map was funded by NASA Grant NAGW-381 2 
and produced in concert by the following: 

Slate of Hewai v i Office of State Planning, 

DLNR Division of Forestry and Wildlife, 

University of HawaTi, Department of Geography 
Geographic Decisions Systems hit* l and Hogan Co. 

Mosaic by Chris Hogan and Leo ikehara. 

Map design by Renee Louis and Eric Yemashita. 
Thanks for contributions by Juho Polo, Steve Sakata, 
Matt McGranaghan, Rob Kiessling, and Harold Gerbiel. 


Lava Flows 

Volcanoes N.P., Big Island 
Change Detection 1 of 6 
From U.H. Geography Dept. 
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Gorse and Koa-Banana Poka 



The image on the left is comprised of 1977 MSS data only. Oor&e and Banana Poka, two exotic 
plants invading the area, are hard to detect in this image. The image to the right is comprised of 
1990 TM data only. Here, Gorse can be seen as a rust color within the yellowish area just below 
the center of the image extending upward toward the keyhole pattern at the top center of the 
image. Koa-Banana Poka is seen as an orange-red. The yellowish areas of the image are the 
high elevation pastures of Mauna Kea. 


Digital data includes 4 scenes from Landset TM (June 1988. 

June 1990. and Oct 1991) registered with 2 scenes from 
Land sat MSS (Oct. 1977) Alt geographic databases of the 
State of Hawaii Geographic Information System are in UTM 
Zone 4, although the Island of Hawai ' i is in UTM Zone 6. This 
image is projected into UTM Zone 4 to comply with the State 
of Hawai ‘ i database design. Base is adapted from Hawai * i Island 
USGS 1:24.000 topographic maps. Old Hawaiian Datum is based 
on Clarke 1 866. MSS processing by Ed Petteys, Chris Hogan, 

Pat Costafes. Mike Buck, and Craig Tasaka. TM processing by 
Chris Hogan. Ron CannareMa, Leo Ikebam, and Roger Imoto. 


This map was funded by NASA Grant NAGW-38 1 2 
and produced In concert by the following: 

State of Hawai ' i Office of State Planning, 

DLNR Division of Forestry and Wildlife, 

University of HawaTi, Department of Geography 
Geographic Decisions Systems Int'L and Hogan Co. 

Mosaic by Chris Hogan and Leo Ikehara. 

Map design by Renee Louis and Eric Yamashita. 
Thanks for contributions by Julio Polo. Steve Sakata, 
Matt McGraneghan, Rob Kiessling, and Harold Garble! 


Gorse and Koa-Banana Poka 
Mauna Kea, Big Island 
Change Detection 2 of 6 
From U.H Geography Dept. 
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Keyhole 



The keyhole pattern in the lower left comer of each image in a Eucalyptus Forest Plantation. The image in the 
upper left depicts 1 977 MSS data of the keyhole. Since the Eucalyptus comprising the rectangular portion of 
the keyhole is not planted it shows up as a different color from the square portion which is planted. The image 
in the upper right depicts the 1990 TM data. Here, the entire keyholo is the same color since the entire region « 
planted. The lower center image depicts both the MSS and the TM data and shows the change between the 
two data as a bright green color. 


Digital data includes 4 scenes from Landsat TM (June 1989. 

June 1990, and Oct 1991) registered with 2 scenes from 
Landsat MSS (Oct 1977). AN geographic databases of the 
State of Hawaii Geographic information System are in UTM 
Zone 4, although the Island of Hawai * i is in UTM Zone 5. This 
image is projected into UTM Zone 4 to comply with the State 
of Hawai i database design. Base is adapted from Hawai' i Island 
USGS 1:24.000 topographic maps. Old Hawaiian Datum is based 
on Clarke 1866. MSS processing by Ed Petteys, Chris Hogan. 

Fat Costales, Mike Buck, and Craig Tasaka. TM processing by 
Chris Hogan, Ron CannereMe, Leo Ikebera, and Roger Imoto. 


This map was funded by NASA Grant NAGW-38 1 2 
and produced in concert by the following: 

State of Hawai ‘ i Office of State Planning, 

DLNR Division of Forestry and Wildlife. 

University of Hawaii, Department of Geography 
Geographic Decisions Systems kvt'L and Hogan Co 

Mosaic by Chris Hogan and Leo ikehara. 

Map design by Renee Louis and Eric Yameshita. 
Thanks for contributions by Julio Polo. Steve Sakata. 
Matt McGraneghan, Rob Kiessling, and Harold Garble! 


Eucalyptus Plantation Keyhole 
Hamakua, Big Island 
Change Detection 3 of 6 
From U.H. Geography Dept 
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Hamakua, Hawai i Island 



B(TM7> GfTMb) »{M5S4) 



TM Landcover Map 


Legend 

Koa. KoaOhia, Ash, or Eucalyptus 
Wet Dense Mixed Forest 
Montane Wet Otuo^Shrh 
Bare Soil 


Legend 

Shadow, Water 

Sugar Cene 

Dry Koa Community 


Legend 

I Montane Wet Exotic Grassland 
CHI Cloud 
■■ Firetree 


| Montane Dry Koa with Banana Poka 
Scale 


10 


10 


Kilometers 


Miles 



1 : 190000 


The top image shows both 1990 TM and 1 977 MSS data, therefore it displays the land use/cover changes 
in this area. The green represents planted fields on both dates. The bright blue represents a change. The 
lower image is the classified map extracted from the TM image with its legend below it. 


Digital data includes 4 scenes from Landset TM (June 1 989, 

June 1 990, and Oct 1991) registered with 2 scenes from 
Landset MSS (Oct. 1977). AH geographic databases of the 
State of Hawaii Geographic Information System are in UTM 
Zone 4, although Ihe Island of Hawai *i is in UTM Zone 5. This 
image is projected Into UTM Zone 4 to comply with the State 
of Hawai ‘ i database design. Base is adapted from Hawai * i Island 
USGS 1:24.000 topographic maps. Old Hawaiian Datum is based 
on Clarice 1 866. MSS processing by Ed Petteys, Chris Hogan. 

Pat Costales. Mike Buck, and Craig Tasaka. TM processing by 
Chris Hogan, Bon CannareMa, Leo Ikebara, and Roger Imoto 


This map was funded by NASA Grant NAGW-38 1 2 
end produced in concert by the following: 

State of Hawai k i Office of State Planning, 

OLNR Division of Forestry and Wildlife, 

University of Rawer i. Department of Geography 
Geographic Decisions Systems Int L and Hogan t 

Mosaic by Chris Hogan and Leo Ikebara 
Map design by Renee Louis and Eric Yamashita. 
Thanks for contributions by Julio Polo, Steve Sakata, 
Matt McGranaghan, Rob KiessKng, and Harold Garble! 


° Change in Cane Cultivation 
Hamakua, Big Island 
Change Detection 4 of 6 
From U.H Geography Dept. 
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Northwest Hawai'i Golf Courses 



This image uses both 1990 TM and 1977 MSS data. A cyan or light blue color represents areas where 
little to no change has occured. The two cyan areas at the top and center of the image have traces of 
green and red which depicts some sort of change. The golf course at the bottom left of the image is 
completely now as it is entirely green and orange. Note the ability to see the golf course design and 
the cloud shadows. The rusty brown area near the cloud shadow is a Buffelgrass-Kiawe mix and the 
yellow brown area in the lower right is Fountain grass. 


Digital data includes 4 scenes from Laridsat TM (June 1 989, This map was funded by NASA Grant NAGW-38 1 1 

June 1990. and Oct 1991) registered with 2 scenes from and produced in concert by the following: 

Landsat MSS (Oct 1977). All geographic databases of the State of Hawai'i Office of State Planning. 

State ol Hawaii Geographic Information System are in IfTM DLNR Division of Forestry and Wildlife. 

Zone 4, although the Wand of Hawai i is in UTM Zone 5. This University of Hawai'i. Department of Geography 

image is projected into UTM Zone 4 to comply with the State Geographic Decisions Systems Ini ' I. and Hogan 

ol Hawai ' i database design. 8ase is adapted from Hawai ' i Island 

USGS 1:24.000 topographic maps. Old Hawaiian Datum is based Mosaic by Chris Hogan and Leo Ikehara. 

on Clarice 1868. MSS processing by Ed Petteys. Chris Hogan, Map design by Renee Louis and Eric Yamashita. 

Pat Costales, Mike Buck, and Craig Tasaka. TM processing by Thanks for contributions by Julio Polo, Steve Sakai 

Chris Hogan. Ron CannareHa. Leo Ikebara, and Roger Imoto Matt McGraneghan. Rob Kiessling, and Harold Gai 


Golf Courses 
Puako, Big Island 
Change Detection 5 of 6 
From U.H Geography Dept. 
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Mauna Kea Area Reserve 


« 




TM Landcover Map 
Scale 


10 0 10 

10 0 

1 : 355777.28 




Legend 

Bare Lava 

Shadow 

Cinder-Ash 

Montane Mesic Exotic Eucalyyptus 
Montane Wet Exotic Grassland 
Montane Dry Open Koe-Mamanc Forest 
Subalpine Dry Open Mamane Forest 
Montane Dry Opm Mamane-Naio Forest 
Subalpine Dry Open Grassland 
Subalpine Mesic Open Shrubiand 


in,-,- aaa Kilometers 
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■ Miles 
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The image at the top left shows only 1977 MSS data. The image at the top right shows only 
1990 TM data. The lower image is the classified map extracted from the TM image with its 
legend to the right of it. 


Digital data includes 4 scenes from Land sat TM (June 1 980* 

June 1 990. and Oct 1991) registered with 2 scenes from 
Landsat MSS {Oct 1977). AH geographic databases of the 
State of Hawaii Geographic Information System are in UTM 
Zone 4 ( although the Island of Hewai * i is in UTM Zone 6. This 
image is projected into UTM Zone 4 to comply with the State 
of Hawai ’ i database design. Base is adapted from MewaTi Island 
USGS 1:24.000 topographic maps. Old Hawaiian Datum is based 
on Clarke 1866. MSS processing by Ed Petteys. Chris Hogan. 

Pat Costales. Mike Buck, and Craig Tasaka. TM processing by 
Chris Hogan. Ron CannareHa. Leo Ikebera, and Roger Imoto 


This map was funded by NASA Grant NAGW-38 12 
and produced in concert by the following: 

State of Hawaii Office of State Planning. 

DLNR Division of Forestry and Wildlife. 

University of Hawai* i, Department of Geography 
Geographic Decisions Systems Int i and Hogan l 

Mosaic by Chris Hogan and Leo Ikeharn. 

Map design by Renee Louis and Eric Yemashita. 
Thanks for contributions by Julio Polo, Steve Sakata, 
Matt McGranaghan. Rob Kiessling, and Harold Garbiel 


Mauna Kea Area Reserve 
Mauna Kea, Big Island 
Change Detection 6 of 6 
From U.H. Geography Dept 
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List of Band Combinations from Landsat MSS and MultidateTM Mosaic 


The MSS mosaic from 1977 and multidate TM mosaic (6/1989, 6/90, 10/91) were 
registered together using DLG's in the State of Hawaii Database. The multidate 
mosaic is projected into UTM Zone 4, Clarke 1866, Old Hawaiian Datum, to comply 
with the state database design. Different spatial, spectral, and temporal resolutions 
are used to reveal change in landscape. 

Bands: TM MSS Viewer Layers: TM MSS 

7-5-4-3-2-I 4-3-2-1 10-9-8-7-6-5 4-3-2-1 


Layer/Band: 

Layer Combinations: 10-2-1 
Band Combinations: 7-2-1 

Layer Combinations: 9-8-4 
Bands Combinations: 5-4-4 

Layer Combinations: 8-4-1 
Band Combinations: 4-4-1 


Layer Combinations: 10-9-4 
Band Combinations: 7-5-4 


Layer Combinations: 10-4-2 
Band Combinations: 7-4-2 

Layer Combinations: 10-9-2 
Band Combinations: 7-5-2 

Layer Combinations: 8-9-3 
Band Combinations: 4-5-3 


Types of Vegetation / Land Cover Change 
Increase in Density of Gorse , Eucalyptus Keyhole 


Mauna Loa Gardens-Mac nut; Hamakua sugar 


Laupahoehoe Open Forests - Banana Poka, Gulches, 
urban growth, Kupaianaha flow 1986-1992 and 
Kalapana, and cloud shadow 

Cloud - opaque, repeat bum area Puna (1977, 1983, 
June 1991); Kupaianaha flow 1986-1992 and 
Kalapana; sugar cane planted both dates, 1984 Lava 
flow from Mauna Loa, 1977 Lava Flow from 
Kilauea, dense biomass increase of Gorse, change in 
wet pasture grasses, and golf coarse and biomass 
increase of Kiawe forest near Mauna Lani Hotel 

Biomass increase in Coastal Kiawe Forest, and Open 
Kiawe Forest - Buffelgrass 

Moisture and biomass increase in Fountain Grass 


Separation of different forest types (i.e.. Myrica faya 
from Ohia, Ohia-Koa from Gorse, with Ohia-Koa 
appearing rusty brown, M. faya and Gorse appearing 
bright red; biomass change in HVNP forest west of 
Kipuka Ki, and separation of vegetation-lava from 
vegetation-soils. 
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The Application of Remote Sensing Data to a GIS Study of 
Land Use, Land Cover, and Vegetation Mapping 
in the State of Hawai'i 



APPENDIX G - Map Composition 

1) Map Composition Process ERDAS to ARC/ Info 

2) Sample of Image Mask Process in ERDAS/Imagine 

3) Large Format Negative Production 



dirdoc.isl 

DESCRIPTION OF FILES UNDER 
/VG03/nasa/finmaps/aml$ 

The two main AMLs in this directory: 

finalmap.aml This AML creates a map of the final supervised classification, 
finalmapraw.aml This AML creates a map of the raw Landsat image. 

The other files under this directory include: 

contour.aml This AML draws the contour lines and is called by finalmap.aml 
and finalmapraw.aml. 

geo.proj Geographic projection file used to draw Lat/Long grid lines. 

utm4.proj UTM projection file used to draw Lat/Long grid lines. 

imgcolor.aml This AML was developed by Royce Jones. It creates shadeset, 

c.map file, and a key file to match a specified grid. A name of a 
grid must be specified in order to run AML. 

DESCRIPTION OF FILES UNDER 
/VG03/nasa/finmaps/covers 

oceanshd This coverage is used to shade the ocean, 

pathrow This is the Satellite Path/Row Index coverage, 

placenames Annotation coverage for placenames. 

DESCRIPTION OF FILES UNDER 
/VG03/nasa/finmaps/grids 

finalgrid.shd Shadeset for final supervised classification grid. 

finalgrid.key Keyshade for final supervised classification grid. 

finalgrid This is the final supervised classifcation grid that was created with 

ARC'S imagegrid command. From this grid, the finalgridclip is 
created to clip out the data in the water. 
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finalgridclip 


bicstgrid 

hawcst 


contour.lkey 
finalinfol .txt 
finalinfo2.txt 
gps.mkey 
rawinfol .txt 
rawinfo2.txt 
rdmang.lkey 
source.txt 


Final supervised classification grid. To view grid in arcplot: 

gridpaint /vg03/nasa/finmaps/grids/finalgridclip # identity ~ 

nowrap /vg03/nasa/finmaps/grids/finalgrid.cnnap 

This grid was created by Royce in order to clip out data in the 

water. 

This is a copy of the Hawaii coastline coverage in 
/vg02/prime/basemap/final/@coast. Royce made a copy of it, 
and defined its coordinate system. 

DESCRIPTION OF FILES UNDER 
/VG03/nasa/finmap$/keyfiles 

Contour keyline file. 

Map information for final supervised classification map. 

Map information for final supervised classification map. 

GPS marker file 

Map information for Raw Landsat map 
Map information for Raw Landsat map 

Major Roads, and Publicly Owned Protected Areas keyline file. 
Map source file. Map source file. 


As of December 27, 1995 the final supervised classification is stored under: 
/vg03/nasa/finmaps/grids/ and it is called finalgridclip. 

As of December 27, 1995 the raw Landsat image is stored under: 

/vg03/nasa/leo/makemaps and it is called hawaiimosaic.lan. 



IMAGE MASKING 

Example 1: 

Water mask - make a copy unsupervised classification image file, select water 
classes and recode to 0, make all other land cover equal to 1. 

Raw image input 1 X Cl water Mask = new raw image for classification. 

Reclassify with water mask to classify only land. 

Example 2: 

Water in IR 7-5 or 4 subset into single channel and recode water-land boundary. 
Separators 0-1; Image Interp>Util>Operators 

Input original image. Input mask X = new raw image for classification. 

Mask a land cover for classification: 

Forest Vis3 or NIR4 same procedures as above. 

Arc Info Polygrid with / without off shore island and new land buffer. 


Appendix C #2 



Island of Hawaii - Landsat Classification Project 
Notes - Rod Low - 2/8/96 

Production of 4" x 5" Film Negatives 

1. Received graphics files from Leo's program. 

2. Used Arc/Info's Postscript Command to produce Encapsulated Postscript 

files . 

3. Used Alchemy version 1.7 to convert EPS to TIFF format 

a. Parameters for Raw Mosaic Map: 

alchemy /gis5/nasa/hawraw. eps /gis5/nasa/hawraw. tif -tl -Zm 2 
-Zi 29 40 -Zo 4i 5i -Zd 3400 3400 -Z+ 

TIF file characteristics: 

Width x Height: 12325 x 17000 

Number of Colours: True Colour (24 bits) 

Dots per inch: 3400 x 3400 

Image size (inches) : 3.62 x 5.00 

LZW Cbmpressed file size: 111 Megabytes 

b. Parameters for Classification Map: 

alchemy /gis5/nasa/hawaii . eps /gis5/nasa/hawaii . tif -tl -Zm 2 
-Zi 40 34 -Zo 5i 4i -Zd 4000 4000 -Z+ 

TIF file characteristics: 

Width x Height: 18824 x 16000 

Number of Colours: True Colour (24 bits) 

Dots per inch: 4000 x 4000 

Image size (inches) : 4.71 x 4.00 

LZW Compressed file size: 38 Megabytes 

4. Negatives created by Snapshot using Adobe Photoshop 

Vectorization of Classified Image (in Arc/Info) 

1. Erdas IMG file converted to Arc/Info Grid using IMAGEGRID command. 

Resulting grid named HAW_SCL. 

2. Coastline vector coverage buffered 1 kilometer. 

3. Buffered coastline converted to Grid using POLYGRID command. 

Resulting grid named COASTBUF . Grid cells out to 1 km from the 
coast received a value of 1, those outside the 1 km buffer 
were given a value of NODATA. 

4. Each quad was processed in the GRID module as follows: 

a. For an area a little larger than each quad, the command 

SETWINDOW * HAW_SCL was given. 

b. QUAD = HAW_SCL * COASTBUF 

c. The GRIDPOLY command was used with a fuzzy tolerance of 1 meter. 

d. Each vectorized coverage was clipped to the boundaries of the quad. 

The item GRID-CODE in the PAT (Polygon Attribute Table) contains 
the value of the original Erdas Image's Class Code. 

Still To do: 

1. Create more generalized coverages. The proposed method is: 

a. Add a numeric item to each quad's PAT table that will signify 
each unique color code from the Erdas Image. 

b. Dissolve each original classified quad coverage using this new item. 

2. Create Lookup Tables with class description text and color symbols. 
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APPENDIX H - Project Products 


The following Products are either too large or on non-paper media to be bound in 
this report. 

1) E-size Color Plot of Big Island Raw Data Mosaic 

- draft plot 

- final on large format 4X5 negative and CD 

2) E-size Color Plot of Big Island Classified Mosaic 

- draft plot 

- final on large format 4X5 negative and CD 

3) Raw Data Final Mosaic Map 

- on large format 4X4 Negative and CD 

4) Vegetation and Other Land Cover Landsat TM Classification Map 

- on large format 4X5 negative 

5) Sample 8X10 size Color Plots of Change Detection 

6) Project CD containing final images and map products 

7) A Pilot Land Cover Classification System for Remote Sensing Use in 
Subtropical Hawaii, based on Landsat TM data. 



